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REPORT OF THE JOINT CSIRO/NASA STUDY 


OPTICAL PROPERTIES ON SOUTHERN HEMISPHERE AEROSOLS 
1. BACKGROUND 

The global measurement of wind is an important aim for meteorolo- 
gists, particularly for numerical weather prediction. Although global winds 
derived from pressure observations using the geos trophic relationship have been 
popular up until the present, it has recently become increasingly clear that 
direct wind data are also extremely effective for use in numerical weather 
prediction. This is particularly true of tropical and southern hemisphere 
regions where there is a paucity of conventional wind and pressure data. 

The realization of global wind measurements from satellites has 
become feasible recently with technological research and development in the 
field of coherent Doppler lidar. Coherent lidar, using visible or infrared 
wavelengths, can provide winds along a cloudless line of sight by measuring the 
frequency change of the laser radiation when backscattered from aerosol 
particles moving with the wind flow. 

The design of a suitable lidar, and even the choice of a suitable 
wavelength, depends on a number of factors such as power and weight require- 
ments and eye safety. The design also depends on a knowledge of the abundance 
of aerosol particles in the atmosphere and their geographical variation 
throughout the troposphere. A global Doppler wind lidar system would obviously 
be of most use in those remote regions of the globe which are poorly served by 
conventional observations. Many of these remote regions are, by definition, 
far from centers of population and therefore have a relatively low level of 
atmospheric aerosol pollution. It is thus imperative for design considerations 
to know the abundance of aerosols, as well as their physical and chemical 
properties, in remote regions of the globe. 

Recently, international programs such as the Global Backscatter 
Experiment (GLOBE) (sponsored by NASA) and the Laser Atmospheric Wind Sounder 
(LAWS) project (also sponsored by NASA) have been set up both to measure the 
global aerosol backscatter and to design a suitable satellite instrument. The 


GLOBE program will operate in the 1988/89 timeframe. 


A preliminary conceptual design for a space- based Doppler lidar was 
undertaken by R. M. Huf faker et al, NOM TM ERL WPL-37 (1978) and NOAA TM ERL 
WPL-63 (1980). The system was named WINDSAT. The apparent feasibility of 
measuring winds from space, based on this study, together with the necessity of 
obtaining aerosol data in "clean" areas of the global atmosphere, led to the 
initiation of collaboration in 1984 between CSIRO Division of Atmospheric 
Research, NASA, and Coherent Technologies, Inc. , a company situated in Boulder, 
Colorado. 


2. SUMMARY 

This study was made in support of the LAWS and GLOBE programs, which 
aim to design a suitable Doppler lidar system for measuring global winds from a 
satellite. Observations were taken from 5°S to 45 °S along and off the E and SE 
Australian coast, thus obtaining representative samples over a large latitude 
range. Observations were made, between 0 and 6 km altitude, of aerosol 
physical and chemical properties in situ from the CSIRO F-27 aircraft; of lidar 
backscatter coefficients at 10.6 Jim wavelength from the F-27 aircraft; of lidar 
backscatter profiles at 0.694 Jim at Sale, SE Australia; and of lidar 
backscatter profiles at 0. 532 Jim at Cowley Beach, NE Australia. 

Both calculations and observations in the free troposphere gave a 
backscatter coefficient of 1-2 x 10” 11 m' 1 sr" 1 at 10.6 ji m, although the 
accuracies of the instruments were marginal at this level. Equivalent figures 
were 2-8 x 10” 9 m" 1 sr" 1 (aerosol) and 9 x 10" 9 to 2 x 10" 8 m' 1 sr" 1 (lidar) at 
0.694 Jim wavelength at Sale; and 3-7 x 10* 9 m" 1 sr” 1 (aerosol) and 10 ” 8 - 10” 7 
m” 1 sr" 1 (lidar) at 0.532 Jim wavelength at Cowley Beach. The measured 
backscatter coefficients at 0.694 Jim and 0.532 Jim were consistently higher than 
the values calculated from aerosol size distributions by factors of typically 2 
to 10. 


It is recommended that much more experimental data be taken and that 
all instrumentation be upgraded, at least to some extent, before further 
measurements sire made, in order to infer the levels of backscatter in the free 
troposphere with sufficient accuracy and sensitivity. 
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3. JO TNT CSTRO/NASA PROGRAM 

The desirability for the development of a viable space program in 
Australia has been recognized recently. In 1984 CSIRO identified participation 
in the development of the Australian space industry as an important strategic 
goal for the Organization, and formed the CSIRO Office for Space Science 
Applications, or COSSA. COSSA has progressed from that beginning to an active 
agency, promoting space research and development between Industry, Government 
Research Laboratories, and Universities. This GLOBE /LAWS experiment was 
sponsored by COSSA and supported by CSIRO, Division of Atmospheric Research. 

The joint field program had four principal aims: 

(i) to make direct measurements of 10.6 fJtm backscatter in the southern 

hemisphere over a wide enough latitude range to be representative and 
during conditions as close as possible to "background". 

Specifically, this meant avoiding periods when dust transport, cloud, 
or other short period enhancements would be likely. 

(ii) to attempt to identify the aerosol component responsible for 10.6 fi m 
backscatter. This was seen as a two part process: 

(a) by determining the aerosol size distribution and using Mie 
theory to calculate backscatter as a function of particle size 
thereby identifying the "active" size range; 

(b) obtain (qualitative) composition data by morphological 

identification of individual particles in the "active" size 
range. 

(iii) to directly measure the ratio of backscatter at 10.6 fxn and 0.693 /im 
and also the ratio of backscatter at 10.6 /im and 0.53 /im to derive 
empirical conversion factors for backscatter at these wavelengths. 
Also, to investigate the aerosol size distribution over the size 
range responsible for this scattering and whether such a conversion 
procedure is in fact valid. 

(iv) to obtain size distribution and species data to add to the very 
limited data set on southern hemisphere free troposphere aerosol. 
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The selected area of operation, near 150 °E and ranging in latitude 
from 4°S to 42 °S was chosen to include conditions ranging from tropical to the 
Southern Ocean. 

The program was planned as a cooperative venture bringing together 
five organizations in two countries. Participants in the field program were: 

- CSIRO Division of Atmospheric Research - airborne microphysics /ground- 
based ruby lidar 

- NASA Marshall Space Flight Center - airborne 10.6 /im CW lidar 

- Curtin University (formerly WAIT) and DSIR - ground-based doubled YAG 
lidar 

Analyses were performed by the individual investigators, with overall 
collation of the data analyses and the analysis of the NASA CW lidar data the 
responsibility of Coherent Technologies, Inc. , Boulder, Colorado, who also 
coordinated the experiment. 

CSIRO, DAR conducted two main experiments. The first involved 
obtaining vertical profiles of aerosol backscatter from the ground and at a 
visible wavelength of 0.694 urn using the CSIRO lidar, and at a site near the 
south-eastern coast of Australia. The second involved flying instruments in 
the CSIRO F-27 aircraft to measure aerosol sizes, numbers and chemical 
properties, at different altitudes. The aircraft could obtain data from 
selected altitudes from the ground up to a maximum of about 6 km. 

NASA participation took the form of a coherent continuous-wave CO 2 
lidar, developed at Marshall Space Flight Center, Huntsville, Alabama. The 
lidar was mounted in the F-27 and looked out sideways. It was able to measure 
the aerosol backscatter coefficient directly at a wavelength of 10. 6 /im and at 
a nominal distance from the aircraft of 10 m. The lidar measured the aerosol 
backscatter coefficient continuously during aircraft experiments. 

The West Australian Institute of Technology, (WAIT) , now Curtin 
University), Department of Applied Physics, also collaborated substantially in 
the experiment. Their participation was part of a program studying lidar 
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backs catter for the Department of Defense, Electronics Research Laboratory. 
They measured backscatter profiles at 0. 532 fim at a tropical location o in 
Northern Queensland simultaneously with the present experiment. 

The Department of Defense, Electronics Research Laboratory also 
contributed financial support, to CSIRO, DAR for the aircraft operations around 
Cairns. 


Coherent Technologies, Inc. , Boulder, Colorado coordinated the 
experiment and also the data analysis. They undertook analysis of the NASA CO 2 
lidar data and also provided facilities for CSIRO to analyze much of the 
aerosol data, and supported one CSIRO scientist for a period during the 
analysis phase. 

A formal agreement based on the proposal was set up between CSIRO, 
COSSA, and the NASA Office of International Affairs, Washington, DC. The 
formal lines of communication and responsibilities of the participating 
agencies are shown schematically in Figure 1. 
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Figure 1. Formal Structure of the Study. 
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4. THE EXPERIMENT 

4.1 INSTRUMENTATION OUTLINE 

4.1.1 A IRB« 

Airborne instrumentation was carried in the CSIRO Space Science and 
Applications Branch (COSSA) Fokker F-27 aircraft. This is a twin turbo-prop 
with a typical cruise speed of 180 knots and ceiling altitude of 6.5 km. 
Support instrumentation included doppler radar and INS for determination of 
wind and location, pitot-static for true air speed and altitude, reverse flow 
temperature sensor, and dew-point hygrometer. 

The main experiment included the NASA-Marshall Space Flight Center's 
10. 6 fim CW lidar system operating through an 18-cm diameter germanium window. 
To conform to the aircraft operational requirements, the lidar orientation was 
perpendicular to its usual disposition requiring some modification to the final 
optics train. 

Aerosol microphysics instrumentation included some real-time devices 
and some particle collectors. A modified GE condensation nucleus counter was 
used for total particle concentrations and in conjunction with a diffusion 
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battery to give size data for radii less than 0. 1 /im. Concentrations of 
particles with radii from 0.05 to 1.5 fim were determined with a PMS ASASP-X 
particle spectrometer. Both the GE and AS ASP were internally mounted and fed 
from the aircraft isokinetic sampling line. A modified CLIMET 208 size 
spectrometer was externally pod-mounted. The normal radius sensing range is 
0. 18 to 10 fim, however an engineering problem arose in the operation of this 
instrument and data from it has not been included in size distribution 
determinations at this point. Particles were collected for later analysis 
using two impaction systems, an external 7-mm wide, free- stream impactor that 
allows up to 5 samples per flight and a 1-rnm jet sampler with morphological 
identification where possible. Because of the problems with the CLIMET 208 
spectrometer, large particle concentrations were determined on representative 
impaction samples despite the intensive manual effort this involves. 

4.1.2 SU RFACE 

The CSIRO ruby lidar was an incoherent lidar transmitting at a 
wavelength of 0.694 fim, with a pulse energy of typically 500 mJ, pulse length 
of 60 ns, and maximum pulse repetition frequency of 1 Hz. The receiving mirror 
was Cassegrain, 30 cm in diameter and aperture of from 1 to 10 milliradians. 
The detector was an RCA photomultiplier. 

The new CSIRO tunable incoherent CO 2 infrared lidar was also 
operated, but due to several engineering problems its range was considered too 
limited to be of any use. 

The Curtin University Nd: YAG incoherent lidar operated at a 
frequency-doubled wavelength of 0.532 /im, pulse energy of 0.2 mJ, pulse length 
of 5 ns, and pulse repetition frequency of 1 Hz. The receiver was of Newtonian 
design and 25.4 cm aperture. 

4.2 WEATHER CONDITIONS 

During the experiment period, weather conditions were typical for 
late autumn-early winter. In southern Australia there was a slow progression 
of anti-cyclonic systems followed by a rapid progression of a series of cold 
fronts across Victoria, after which deep south or southwesterly streams 
penetrate into continental southeastern Australia. Measurements at the Cape 
Grim baseline station in northwest Tasmania show that these deep southwest 
streams usually produce clean or "baseline" conditions for both gases and 


8 



condensation nuclei. The frequency of baseline conditions varies from year to 
year, but is typically 30 to 40% (depending on definitions). Weather 
conditions for the northern section of the experiment were also quite typical 
for the season with low level southeasterly trades giving way to the upper 
level weatherly zonal flow above the inversion. 

4.3 LOCATIONS AND FLIGHT PLANS 

This study was planned around two principal sites, each with a 
ground-based lidar and a series of "excursion" and "transit" flights. The two 
fixed sites were at Sale (38°06 / S, 147°04 / E) in eastern Victoria where the 
CSIRO ruby lidar was located, and Cowley Beach (17°41 / S, 146°07 / E) in northern 
Queensland where the Curtin University doubled YAG lidar was sited for the 
experiment. In both of these locations, aircraft flights were planned as a 
series of straight and level legs, each of approximately 20 to 30 minutes 
duration centered on the lidar site. Altitudes were stepped between the 
boundary layer and the aircraft ceiling (usually between four and six levels). 
Flights normally finished with a spiral descent sounding. Excursion flights 
south from Sale to the Cape Grim (40°41 / S , 144° , 41 / E) area, and north from 
Cairns to Rabaul (4°12 / S, 152° 12 7 E) in New Guinea were included to extend the 
geographic range of measurements. Both the excursion flights and transit 
flights used stepped levels in the free troposphere to obtain representative 
samples (see Figure 2). 


9 










5 . 


DATA ANALYSIS - SUMMARY OF MEASURED DATA 

Details of the data analyses are given in the various appendices. 
Only the final analyzed data are presented here. 

5. 1 AEROSOL 

Details of the data analysis and the results are given in Appendix A. 
This section surrmarizes the most significant results. Profiles of 
concentration for particles with radii larger than selected values plotted as a 
function of altitude showed several important features. Concentrations of 
optically important particles were relatively uniform or declined slowly with 
increasing altitude in the boundary layer, but above the inversion, decreased 
steadily to an altitude of about 2.5 or 3 km above which concentrations were 
relatively constant. The larger the particle size the more rapidly 
concentrations fell off with altitude above the mixed layer. In Sale, the 
lower inversion meant that free troposphere conditions were approached at a 
lower altitude than at Cairns and also the concentration of large particles , 
for example, those with radii larger than 0.2 /on, was somewhat lower there 
also. Total particle concentrations indeed confirmed a very clean air mass 
during the period at Sale when large particle concentrations were low, and more 
typical conditions during the observations near Cairns. 

Particle size distributions were obtained as means over selected 
altitude bands and geographic regions. The measured distributions followed the 
well known multi-modal form but typically only exhibited two main modes. In 
the lower layers (below 2.3 km) the individual morphology of the large 
particles (r > approx. 1 fim) was predominantly a sea-salt type; however, the 
recrystallization did not produce as regular cubes as is usual. Approximately 
12% of particles for these altitudes were irregularly shaped and could not be 
classed in this way. Low level winds were either maritime or continentally 
modified maritime during the experiment. 

Most free troposphere particles from both main sites, when examined 
by electron microscopy, had a morphology typical of sulfuric acid or lightly 
aranoniated sulfate. Only 24% of the free- troposphere particles examined didn't 
have this morphology. 

Day to day variations in particle concentration are an important and 
regular feature of the atmospheric aerosol. During this flight series, the 
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most obvious variations in the free troposphere occurred at Sale in a 
post-frontal, deep southwesterly flow where the concentrations of large 
particles were about one- tenth of the average (for the whole period) . This was 
followed by an increase of about two orders of magnitude (at some altitudes) 
with the reversion to a continental airstream. 

5.2 MIE ANALYSIS 

Full details of calculation and analysis are given in Appendix A. 
Computation of backscatter coefficients from the aerosol size distribution data 
serves two purposes. First, by calculating and plotting differential 
backscatter distributions, e.g. as log (dlog b/dlog r) as a function of log r 
the 'active' aerosol, that is the fraction contributing significantly to the 
backscatter at a given wavelength is readily identified. Second, the 
calculated value of backscatter over the whole distribution gives a means for 
comparing experimentally determined backscatter coefficient values and 
experimentally determined aerosol concentration values. 

Backscatter values were calculated for a series of specific bimodal 
size distribution models corresponding with the chosen altitude bands, and at 
wavelengths of 0. 53 fim, 0. 693 fan and 10.. 6 ftta. Refractive indices were selected 
to encompass a likely range of compositions. The first assumption is that the 
aerosol mode centered at 0.06 fim radius (the ubiquitous "accumulation' mode) is 
sulfate, however the degree of neutralization is unknown and a wide range of 
cat- ion composition and hydration combinations is possible. As a 
simplification this was assumed to be ammonium bi-sulfate. More evidence is 
available on the composition of the large mode particles. In the boundary 
layer and just above the inversion (z = 2.3 km) irrpaction samples indicated a 
high proportion of sea-salt particles. In contrast, the upper troposphere 
samples showed predominately sulfate and only partially ammonia ted. For the 
lower levels two model compositions were used for the large mode component; 
'maritime' or water. Since the refractive index of sea-salt is quite strongly 
dependent on the actual relative humidity, particularly for RH > 60%, this 
approach should span the likely range of refractive indices. 

Separation of the two aerosol modes for the calculations demonstrated 
clearly that at 10. 6 fi m the backscatter was dominated by large mode particles 
for all altitude bands and regardless of the refractive index model used. In 
the free troposphere the contribution fell to 10% at radii of 0.3 fi m and 2 fm, 
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peaking at around 0.8 /im radius. In contrast, at shorter wavelengths, 
particles in the smaller 'sulfate' or 'accumulation' mode made an equivalent or 
dominant contribution to the backs catter compared with the large mode. It 
follows that backscatter at these widely different wavelengths may be dominated 
by particles with different sources, sinks and histories. In the boundary 
layer at all the wavelengths considered, backscatter was dominated by the large 
particle mode with radii from around 0. 5 (ten and up to 10 or 20 ft m being 
important in the integrated value. 

5.3 AIRBORNE CW C0-> LIDAR 

Full details of the analysis and results of the CO 2 lidar aerosol 
backscatter measurements in Australia are given in Appendix B. Plots of 
measured aerosol backscatter with altitude are presented for - - 0 to 6.5 km 
altitude, 5° S to 42° S latitude, and 145° to 153° E longitude. Measurements 
with the NASA-MSFC CW CO 2 lidar in the free troposphere gave aerosol 
backscatter coefficient values of 1-2 x 10 _1 1 m -1 sr _1 at 10.6 fm. The aerosol 
backscatter in the free troposphere at 10.6 /im did not change significantly 
from 5° S (Rabaul) to 42° S (Tasmania). The accuracy of the lidar was marginal 
at this level. The minimum sensitivity of the lidar was found to be 
approximately 2 x 10 -11 m“ 1 sr~ 1 . The lidar data was taken in two modes of 
operation; a volume scattering mode and a single particle counting mode. Only 
the volume scattering mode data was analyzed. For the CW CO 2 coherent lidar at 
10.6 fim, particles around 1 fi m were the most important due to the rapid falloff 
in number density above 1 fan. 

It is recommended that more experimental data be taken in Australia to 
measure the data base in this very "clean" southern hemisphere environment. 
The CW CO 2 lidar needs to have increased sensitivity to levels of 5 x 10” 12 
m _1 sr" 1 . An increase in laser power seems the most effective way to increase 
the CW lidar aerosol backscatter sensitivity. 


5.4 SURFACE LIDAR 

5.4.1 CS IRQ /PAR RUBY LIDAR AT SALE 

Full details of the analysis and results are given in Appendix C. 
The final data are presented as height profiles of aerosol backscatter 
coefficient representing different time and height averages. The main features 
in the backscatter were the generally high, but variable values below the mixed 
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layer inversion, followed by the rapid decrease by several orders of magnitude 
to a minimum value which lay at altitudes between 1.2 km and about 4 km. On 
one day when the air above the inversion was essentially of continental origin, 
a very sharp minimum at 1.2 km was followed by a recovery in aerosol 
hacks catter which then remained fairly constant with altitude. 

The backs catter coefficient in the mixed boundary layer varied by 
mere than an order of magnitude from day to day, with high values immediately 
following an incursion of racist southwesterly air, and lower values in air 
which had a mere continental origin. 

5.4.2 CURTIN UNIVERSITY Nd: YAG LIDAR AT COWLEY BEACH 

Full details of the analysis and results are given in Appendix D. 
Lidar data are presented as height profiles of aerosol volume backscatter 
function (at 532 nm) for the 9th and 10th of June 1986, the two days on which 
useful comparisons could be made between the ground-based lidar and the 
airborne particle samplers. Both profiles were averages of 50 lidar firings 
and showed enhanced aerosol scattering in the mixing layer, with maximum values 
of 5 x 10~ 7 m -1 sr -1 to 2 x 10 ” 6 m -1 sr" 1 . Above the temperature inversion, 
aerosol backscatter ing , as measured by the lidar, decreased steadily with 
height to values of 10" 8 m -1 sr -1 or less. Aerosol profiles were normalized to 
molecular scattering values at around 4.5 km altitude, so that any significant 
aerosol scattering at these levels meant that the lidar data underestimated the 
aerosol scattering at all other levels. Because aerosol scattering was weak 
around these levels, any underestimate is expected to be minimal. 

As explained in Appendix D, a single measured value of the aerosol 
backs catter-to-extinction ratio was used in the calculation of the 
transmittance correction during data analysis. If this value was not constant 
with height or differs from the value used, the profiles of backscatter would 
have to be scaled accordingly. 

5.5 UNCERTAINTIES. DATA RELIABILITY 
5.5.1 AEROSOL 

From the outset it should be appreciated that sampling large aerosol 
particles from aircraft is difficult. Concentrations may be very low, 
dictating large flow rates and/or long sample times to obtain good estimates of 
concentration. Furthermore, there are potential losses or gains in 
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concentration in sampling non-isokinetically and due to the effects of flow 
compression and expansion around the aircraft as well as screening off certain 
areas to particle flow. Plumbing particles along sampling lines with 
decelerator stages, bends, stepped sizes, or joins may involve particle losses 
through turbulence or impaction. These effects will vary with particle size, 
altitude and aircraft speed. In general, quantifying the magnitude of such 
losses if possible, is exceedingly difficult. 

Light-scattering, particle-size spectrometers are subject to 
uncertainties in size determination relating to the accuracy of the calibration 
relation, unknown refractive index of sample particles, conversion from 
calibration particle refractive index to assumed refractive index of sample 
particles, dependence of reponse on particle shape, and finite response width 
of the counter (particularly where the distribution slope is very steep) . 
Snail errors in size can appear as large concentration errors. A reasonable 
estimate of the sizing error if the composition is known is around 5 to 10%, 
but this could be much larger with unknown composition (and if no conversion 
from a PSL calibration size is used) . Other techniques such as impaction 
require a model to convert the impacted, dess ica ted electronmicrographed image 
back to an equivalent ambient particle size. If the irnpactor operates in a 
region where the collection efficiency is less than unity, the uncertainty in 
the collection efficiency, which is usually determined empirically, must also 
be included. All possible efforts were taken to reduce the potential for 
errors in sizing and counting, as an overall estimate sizing error should lie 
in the range of around ± 5% for particles with radii around 0. 1 fJtm to about ± 
20% for radii greater than about one micrometer. 

5.5.2 AIRBORNE CW COs LIDAR 

Techniques were developed during this study effort, and described in 
Appendix B for determining the minimum sensitivity of the CW CO 2 lidar by 
analyzing the noise records of the lidar data. There has been considerable 
effort at NASA-MSFC to develop the CW CO 2 lidar and to develop proper 
calibration techniques. The volume scattering mode is more understood than the 
single-particle counting mode and was used in the data analysis. However, the 
minimum sensitivity of the CW CO 2 lidar was found to be approximately the same 
value as the measured data in the free troposphere. The CW CO 2 lidar needs 
additional laser power to increase the lidar sensitivity to 5 x 10“ 12 m~ 1 sr~ 1 , 
as well as special care in calibration. The single- particle mode needs to be 
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better understood and calibrated for future testing. Improvements can also be 
made in using calibration targets when in flight and by verifying that the 
processing bandpass filter is correctly matched to the true airspeed of the 
aircraft. The Dicke switching technique is subject to true airspeed errors, 
pointing angle errors, and filter tuning errors. Techniques should be 
developed to confirm that the processing bandpass filter setting is correct. 
Quick- look analysis procedures should also be developed to obtain a near 
real-time assessment of the data quality. The data gathered in Australia 
indicated that the 10.6 fi m aerosol backs catter in the free troposphere is of 
the sane level as the noise of the lidar. It is recommended that similar 
flight experiments be again conducted in Australia, but with increased CW CO 2 
lidar sensitivity, improved system calibration, and an on-board quick- look 
analysis procedure. 

5.5.3 SURFACE LIDAR 

The measurement of aerosol backs catter coefficient in the visible 
spectrum and in the free troposphere is difficult in several respects. First, 
in the clean maritime atmosphere of the Southern Hemisphere, the aerosol 
backs catter coefficient in the free troposphere is found to be typically one to 
two orders of magnitude less than the molecular backscatter , particularly at 
the shorter wavelengths. Thus, errors involved in extracting a (usually) 
slowly varying aerosol signal from the molecular signal, and the associated 
normalisation problem, are considerable. Second, the (range) ^ decrease in 
lidar signal implies that signal- to- noise ratios required to extract the 
aerosol signal from the molecular signal are marginal, at best, in the upper 
free troposphere, unless many hundreds of profiles are averaged. Third, offset 
errors due to background radiances (a daytime problem) and due to amplifiers in 
the signal line can cause appreciable, consistent departures from reality in 
the aerosol signal. 

In many respects, the accurate measurements of the free tropospheric 
aerosol backscatter coefficient in clean atmospheres requires techniques 
similar to those used for many years in the stratosphere. In particular, 
photon counting methods would appear to have some advantages over analog 
amplification techniques. 
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6. 


COMPARISON OF MEASURED AND CALCULATED BACKSCATTER COEFFICIENTS 

Backscatter coefficients were calculated as a function of altitude 
for the three wavelengths used in the field experiment from the mean aerosol 
concentration profiles measured with the ASASP and measured (altitude dependent 
model) aerosol particle size distributions. Backscatter coefficients 
calculated for the model size distributions were scaled for different altitudes 
using the ratio of the number of particles larger than a specified radius (0. 06 
/m or 0.25 frn) in the rrean concentration profiles and the corresponding 
concentration in the mean distribution. This gives continuous profiles of 
calculated backscatter with altitude but with breaks at the altitude boundaries 
where the size distribution model changes. The profiles for 10.6 /im, 0.694 /im, 
and 0.532 /on backscatter calculated in this way are given in Figures 3 to 5, 
respectively, with experimentally observed profiles of backscatter. 

The comparison for 10.6 /im backscatter shown in Figure 3 differs from 
the shorter wavelength comparisons with respect to the data sets actually 
compared. For 10.6 /im both aerosol and backscatter data apply to the whole 
experimental period, although the calculated profiles were derived from aerosol 
averages (and selected size distributions) in the three different geographic 
regions. For the shorter wavelength profiles, a mean experimental profile is 
plotted for 0.694 /im which comprises data from four of the seven days that are 
included in the aerosol data, and for 0.532 fsm only two profiles (two days) 
from the seven days with aerosol data in the Cairns region are included (at 
this time). This complicates the comparisons to some degree, however there are 
a number of conclusions that can be drawn. 

Considering first the 10.6 ftm data, it is clear from the distribution 
of data points and their associated errors bars in Figure 3 that the 10.6 /im CW 
lidar was operating just on or in its noise threshold in the free troposphere 
for most of the conditions sampled. With the reduced complement of particle 
measuring equipment , the sane can be said for the 10. 6 /«n backscatter. Despite 
this the agreement between the observed and calculated backscatter at 10.6 /im 
is surprisingly good. The backscatter profile was calculated (as above) from 
the concentration (profile) of particles with radii greater than 0.25 /im and 
also assumed a free troposphere composition of ammonium sulfate and boundary 
layer composition of ammonium sulfate in the accumulation mode and water in the 
large particle mode. In the transition between boundary layer and free 
troposphere, a maritime refractive index was used for the large particle mode. 
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The more rapid fall off in observed backscatter contribution is strongest at 
10.6 fim for particles with radii around 1 fim. Above about 1 km, the 
concentration ■ of these particles falls off with increasing altitude 
significantly faster than the 0.25 fim particles on which the profile was based. 
The peak in backscatter at about 4.5 km in both the observed and calculated 
profiles is an artifact due to the presence of widespread altostratus cloud 
north of New Guinea. 

Perhaps the greatest difficulty with analysis of data from both of 
the shorter wavelength lidars, which were operating in daylight conditions, has 
been establishing the absolute calibration. In both cases the aerosol 
backscatter in the free troposphere was only a snail fraction of the molecular 
backscatter. Because of this and also because of the (range) 2 effect on gain, 
the signal due to aerosol backscatter could not be resolved in the random noise 
at the upper levels. For comparisons at both wavelengths, two calculated 
backscatter profiles are given in Figures 3 and 4 for each location. 
Backscatter coefficients were again determined from the ratio of the mean 
concentration of particles with r > 0. 06 fan (continuous line) and r > 0.25 fan 
(dashed line) determined with the- ASASP and the corresponding concentration in 
the model size distribution used for that altitude with the backscatter 
coefficient calculated for the mean size distribution. Calculations of the 
contribution to backscatter as a function of particle size in Appendix A, 
indicated that for both of these wavelengths the major contribution to 
backscatter in the free troposphere and about 25% in the boundary layer was due 
to particles in the accumulation mode. On this basis, the concentration 
profile for particles with radii greater than 0. 06 fan should give a better 
representation than the 0.25 Jim profile. In comparing the observed and 
calculated 0.694 fan backscatter profiles for Sale, Figure 4, backscatter 
calculated using the larger than 0. 06 fan radius particle profile agrees more 
closely with the observed profile between about one and five kilometers 
altitude, but both give considerably smaller backscatter coefficients than was 
observed directly at around 1 km altitude. However, the direct observation of 
backscatter coefficient at 1 km was nearly an order of magnitude greater than 
on any other days. Above about 4 km there is an increasing divergence between 
the observed and calculated backscatter coefficients. This could arise either 
from background and range corrections in the lidar or could represent an 
increasing loss of some particle fraction with increasing altitude; however, 
for this short wave scatter the principal contributing size fraction is not the 
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roost difficult to measure. To some extent here also, the lidar data was biased 
upwards by two profiles on May 25 where the backs cat ter was consistently higher 

C 

in the free troposphere than on the other days. 

Comparison of the observed and calculated variation of 0.532 ftm 
backscatter coefficient with altitude shows an evidently even greater 
discrepancy than was observed at 0.694 /an, although as explained above, only 
two observed profiles have been included at this time. As well, the calculated 
profiles include a number of days where the concentrations were relatively low 
compared with the days of the directly measured profiles. Despite this, there 
appears to be a systematic difference in the shape of the observed and 
calculated mean profiles above the inversion with the aerosol data indicating a 
rapid drop in backscatter- coefficient above the inversion and a much mere 
gradual drop in the lidar profiles. The rapid drop in the lidar profiles at 
about 4 km indicates a close proximity to the height of normalization to the 
molecular profile. 

In conclusion, it was found that measured visible wavelength 
backscatter coefficients were consistently greater than values calculated from 
aerosol size and concentration measurements but that agreement between 
observation and calculation was better than 10.6 /im. Thus, assuming that the 
shorter wavelength lidars are correctly normalized (and therefore calibrated 
correctly) , they are either "seeing" particles which the aerosol 
instrumentation doesn't detect or there nay be a consistent undersizing of 
particles. 


7. RECOMMENDATIONS 

The performances of all three lidars are marginal for measurements of 
the free troposphere Southern Hemisphere aerosol. It is thus recommended that 
effort is put into improvement of these instruments before further observations 
are made. Specifically: 

LIDAR INSTRUMENTATION 

7.1 For the CW 10. 6 /im lidar, a greater sensitivity, that is, a better 
signal- to- noise ratio for a given aerosol backscatter coefficient is required. 
The calibration of the instrument should include a measure of the real zero. 
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7.2. Lidar experiments in the visible region should be done at night. 
Many profiles integrated over an extended time are required. Photon counting 
methods of signal recovery would be preferable. Some form of independent 
calibration is required. 

7.3 There is some difficulty in tying up the relationship in reconciling 
measured and calculated values of backs catter coefficients, particularly in the 
visible region. Support instrumentation to measure extinction coefficients, 
such as nephelometers , would be desirable. 

AEROSOL INSTRUMENTATION 

7.4 As wide a range of instrumentation as possible is desirable, 
particularly instruments that can be mounted externally and have a minimum 
disturbance on airflow. For calculating 10.6 /an backscatter, particular 
emphasis should be placed on particles in the size range 0. 5 fim to 50 pm. 
Further, the composition of the particles for this size range should also be 
determined. 

7 . 5 Measured aerosol concentrations in the Southern Hemisphere appear to 
cover a wide range of magnitudes , both spatially and temporally. Many more 
measurements are needed before a reasonably accurate assessment of aerosol 
backscatter probabiity distribution can be made. 
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Figure 3 . Measured 10.6yim backscatter coefficients from the NASA CW lidar 
for the whole experiment, cloud-free conditions (points with error bars). 
Calculated values based on number of particles with r>0.25iim are for Sale 
(dash line) , Cairns (solid line) and New Guinea (dot-dash line) . 
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Figure 4 . Mean profile of measured 0.694pm backseat ter coefficient from 
Sale (light continuous line) and two calculated profiles derived using 
averaged size distributions normalized by the concentration of particles 
with r>0.06pm (continuous heavy line) and r>0.25pm (dashed line), from the 
mean of all ASASP data for Sale. 
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Figure 5 . Two observed 0.532pm backscatter coefficient profiles from 
Cairns (dot-dash line) and two calculated profiles using averaged size 
distributions normalized by the concentration of particles with r>0.06pm 
(continuous heavy line) and r>0.25pm (dashed line), from the mean of all 
Cairns area ASASP data. 
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INTRODUCTION 


This report describes results obtained from a series of airborne aerosol- 
particle-microphysics measurements made during a joint CSIRO-NASA field 
experiment in the eastern Australia-Papua New Guinea region, between May 
25 and June 11 1986. These aerosol measurements were made concurrently 
with (and from the same aircraft as) 10.6pm backscatter determinations by 
W. Jones of NASA MSF. During the same period backscatter was also measured 
at two shorter wavelengths, 0.693pm and 0.532pm, using ground-based lidars 
sited at Sale in south-eastern Victoria and Cowley Beach in north 
Queensland respectively. Descriptions of the lidar systems, their methods 
of operation and results are described elsewhere. The altitude range 
covered during the airborne measurements was surface to approximately 6km. 

EXPERIMENT DETAILS 
Instrumentation 

All of the airborne instrumentation was carried in the CSIRO Fokker F27 
research aircraft. This is a high-wing twin turbo-prop, aircraft with a 
nominal ceiling of 6.5km. Environmental data recorded during the flights 
using the microphysics instrumentation were drybulb temperature, using a 
reverse flow thermometer, dew-point temperature using a Bendix cooled 
mirror hygrometer and pitot and static pressures, for air-speed and 
altitude. Position information and wind speeds were derived from INS and 
doppler radar systems operated by the aircraft facility. 

Aerosol microphvsics instrumentation 
Summary 

1. Modified GE condensation nucleus (CN) counter with six stage screen 
type screen diffusion battery. Radius range approximately 0.002pm to 
0.1pm. 

2. PMS ASASP-X active cavity laser single particle size spectrometer, fed 
from aircraft isokinetic sampling line. Radius range nominally 0.05pm to 
1.5pm. 

3. Climet 208 white light single particle size spectrometer, modified and 
mounted in external pod. Two selectable inlet nozzles, one isokinetic the 
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other sub-kinetic. Two radius ranges 0.2pm to 1pm and 0.5pm to 10pm. 

4. External, free-stream impactor. A pneumatically operated five slide 
impactor used for collecting particles with r > apx. 2pm on 7mm by 25mm, 
silicone grease coated slides (mounted on 7 mm diameter support rods). 

5. A 1mm diameter round jet, single stage inertial impactor used for 
collecting particles with r > 0.3pm (apx.) at a nominal flow rate of 4 1 
min -1 . Particles were collected on standard, nitro-cellulose coated, 3mm 
diameter electron microscope grids further strengthened with a thin film 
of carbon. 

Operation of aerosol sizing instrumentation 

Aerosol particle concentrations were determined in-situ as a function of 
particle size using a PMS ASASP-X size spectrometer fed from the aircraft 
isokinetic sampling line. This spectrometer covers a nominal radius range 
of 0.05pm to 1.5pm in four overlapping groups of fifteen size channels, 
however the sampling rate is quite low, typically 1.5 cm^s - *. It is 
usually necessary to integrate for considerable periods to obtain 
meaningful concentrations for particles with radii greater than a few 
tenth micrometre. Normally, level flight legs were about thirty minutes 
duration giving around 2.7 litre samples with the ASASP per altitude 
section per flight. Calibrations with PSL particles were carried out 
before, during and after the flight series and corrections made for 
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expected refractive index based on the work of Garvey and Pinnick A . 
Throughout the flight series the ASASP performed well. One problem that 
appears to be common to this family of instruments is a tendency to count 
excess particles in the most sensitive channel when operating at altitudes 
above about 5km and using the recommended ratio of sample to purge air. 
This happened occasionally in the early flights but in general was avoided 
by operating at a flow ratio of around 1 to 15 compared with the 
recommended 1 to 20. There are some indications of a possible inertial 
loss problem somewhere in the sampling-counting process. As it will be 
shown later concentrations of particles measured with the ASASP for radii 
greater than around 1pm appear to fall too rapidly (with increasing 
radius) . This effect is only evident in the boundary layer where the 
concentration of particles with r > 1pm is large enough to measure (with 
the ASASP) and there is another measurement of concentration independent 
of the isokinetic sampling line. 
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For this flight series the main instrument for sizing and counting 
particles with radii greater than about 0.2pm was to have been a 
(modified) Climet 208 externally (pod) mounted and using either an 
isokinetic or subkinetic (concentrating) inlet nozzle. This instrument 
operates at a nominal flow rate of 7 1 m -1 and can be used to greater than 
10pm radius. Examination of the PHA subsequent to the flight series 
however revealed an incorrect internal setting that rendered data for 
radii larger than about 1pm totally irretrievable and greater than about 
0.5pm likely to be low in concentration. Consequently data from the Climet 
has been excluded entirely from size distribution determinations although 
it is used to show comparitive altitude variation and day to day changes 
in concentrations. 

For smaller particles (typically r < 0.1 pm ) penetrations through a 
screen diffusion battery^ were determined with a modified GE type 
condensation nucleus counter^. These penetrations were inverted using 
Twomey's non-linear iterative scheme^ to derive size distributions in the 
radius range 0.002pm to 0.1pm. These diffusion battery determinations were 
initiated manually during the flight at fairly regular intervals and for 
the remaining time the condensation nucleus counter was used to monitor 
the CN concentration. 

Larger particles were collected with a free-stream (externally mounted) 
remotely operated impaction system on 7mm-wide silicone grease coated 
slides (principally in the boundary layer) and on 3mm diameter electron 
microscope grids using a 1mm diameter jet impactor fed from the aircraft 
isokinetic sampling line at a flow rate of 4 to 5 1 min - *. These 
collection systems were included principally to obtain compositional data 
but were used for size information also because of the failure of the 
Climet spectrometer PHA. Collected particles were sized by first 
photographing slide cross-sections (or grids) in either secondary emission 
or transmission scanning electron microscopy, measuring the particle 
images and converting the dry impacted volume to an equivalent 
preimpaction spherical size. Where appropriate, further correction was 
made to account for ambient humidity using the relations given by Hanel^. 
For sizing sulfate particles in the free troposphere the model of Gras and 
Ayers^ was used. Salt particles were modelled as "blocks"; for samples 
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measured using transmission microscopy the height was determined by 
measuring the length of individual particle shadows obtained from a non- 
reactive coating evaporated at a known angle onto the sample before 
microscopy. The mean ratio of height to width obtained in this way was 
also used for samples measured using secondary emission microscopy. 
Particles that could not be identified morphologically were treated as 
spheres with a diameter equal to the mean of the measured dimensions. 
Collection efficiency for the jet impactor has previously been determined 
experimentally^ and for the slide impactor an empirical determination by 
Wong et al. 8 was used. In both cases only particles for which the 
collection efficiency was greater than 50% were included in the analyses 
(after correcting the concentration for collection efficiency) . 

Locations and flight plans 

This study was planned around two main sites, each with a ground based 
lidar system and included a series of supporting "excursion" and "transit" 
flights. The two sites were Sale (38 o 06'S,147 o 04'E) in eastern Victoria 
where the CSIRO ruby lidar was located, and Cowley Beach 
(17 o 41'S,146°07'E) in north Queensland where the Curtain University 
doubled YAG lidar was located, both lidars being temporarily sited for the 
experiment. In both of these locations aircraft flights were planned as a 
series of straight and level legs, each of approximately thirty minutes 
duration centred on the lidar site. Locations of the main sites, local, 
transit and excursion flights are all given in Fig.l. Altitudes were 
stepped between the boundary layer and the aircraft ceiling (usually 
between four and six levels) . Flights normally included a spiral descent 
sounding. Excursion flights south from Sale to the Cape Grim area (40° 41' 
S, 144° 41* E) and north from Cairns to Rabaul (4°12'S,152°12 , E) in New 
Guinea were included to extend the geographic range of the study. 

Excursion flights and transit flights were made at a series of altitudes 
in the free troposphere to obtain a range of representative data. 

Weather conditions 


During the experiment period weather conditions were fairly typical for 
late autumn early winter. In south-eastern Australia the dominant feature 
was the slow progression of an anti-cyclonic system and rapid progression 
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of a series of cold fronts across Victoria. Deep south-westerly streams 
associated with these fronts usually produce clean or "baseline" 
conditions for both gases and condensation nuclei. The frequency of the 
"baseline" conditions varies from year to year but at Cape Grim in north- 
western Tasmania (see Fig. 1) is typically 30 to 40% (depending on 
definition) . Weather conditions for the northern section were also 
typical for the season with low level south easterly trades giving way to 
the upper level westerly zonal flow above the inversion. 

OBSERVATIONS 


Because aerosol particle concentration (and composition) is a function of 
size, altitude, location and time there is no one simple method to fully 
depict the state of the aerosol. The two procedures followed here will be 
to give representative profiles of the number concentration of particles 
larger than a given size, as a function of altitude for particular 
locations and times and to give averaged size distributions covering 
certain times, geographic locations and altitude ranges. 

Concentration profiles were usually determined during a spiral descent but 
also occasionally during a normal ascent. Figures 2 to 5 show a series of 
representative profiles combining output from the GE CN counter and 
cumulative concentrations from the ASASP and Climet size spectrometers. 
Only two lower radius bounds from the ASASP (channels 2+ and 21+ giving 
the concentration of particles with radii greater than apx. 0.06pm and 
0.13pm) and two from the Climet (channels 2+ and 6+ giving the 
concentration of particles with radii greater than apx. 0.2pm and 0.5pm) 
were used. In the configuration used the GE counter gives the 
concentration of particles with radii greater than about 0.002pm. 

The first of these profiles, Fig. 2, a descent into Sale on May 29 1986 
illustrates an important feature observed during the period in Sale. 
Following the passage of a cold front a deep south-westerly air stream 
bringing air from over the southern ocean crossed Victoria, see Fig. 6. 
Particle concentrations in this airstream were exceedingly low . In Fig. 2 
two altitude regions, 0.7km-1.7km and about 3km-5km show total particle 
concentrations below 100 cm -3 and large particle concentrations are also 
severely reduced in the upper altitude band. With the normal descent rate 
of 300 m min -1 several of the ASASP records (one minute sample rate) 
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contain no counts for radii greater than ,13pm, the Climet also shows 
marked depression of concentrations of particles with radii greater than 
0.2pm. Below the inversion (at 1.8 km) large particle concentrations are 
relatively constant with altitude but there is a strong decrease above the 
inversion. Cloud was thin strato-cumulus (apx. 300m thick) with tops to 
around 1.7km. The following day the eastwards progression of the 
anticyclonic system centred in the Great Australian Bight brought upper 
level air that had passed over central Australia to south-eastern Victoria 
and northern Tasmania, see Fig. 7. Particle profiles illustrated in Fig. 3 a 
descent off Cape Grim on May 30 show that there is still extensive 
scavenging of small particles evident below the inversion with total 
particle concentrations approaching a few per cubic centimetre and 
essentially no particles with radii less than 0.06pm but above the 
inversion there is now a marked increase in the concentration of large 
particles with the concentration of particles with r>0.13pm at 3.5km 
greater than observed anywhere else in the study except one descent in New 
Guinea. Comparison of Figures 3 and 2 shows in fact that in the region 
3.5km-4km there is a change in concentration of approximately two orders 
of magnitude in the number of particles with radii greater than 0.2pm (and 
greater than 0.13pm) over the period of one day in the same geographic 
region. 

Profiles for the descent into Rabaul on May 4 are shown in Fig. 4. Although 
there was considerable cloud in the area (cumulus with tops > 6km and 
altostratus at 4.5km) the major part of the descent was cloud-free with 
one cloud layer at 3.3km and some rain at about 1km. As this descent was 
early afternoon, mixing was reasonably well developed. By comparison a 
subsequent morning ascent (see Fig. 12) showed concentrations at 2km about 
two orders of magnitude lower than during this descent. The CN profile for 
this descent is interestingly flat with altitude outside the immediate 
boundary layer where there was a considerable amount of local smoke. 

June 10 was selected to represent soundings taken at Cowley Beach in north 
Queensland and indeed is quite representative of the whole period. The 500 
hPa stream analysis given in Fig. 8 shows that the immediate origin for 
air arriving at Cowley Beach on June 10 was over central Australia. 
Following the air stream back several days shows the origin to be to the 
north west of Australia. Large particles show a levelling out of 
concentration above about 2.3km the altitude of the inversion. During the 
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northern New Guinea show up as a large particle enhancement at about 4.5 
km altitude in the New Guinea profiles. 

Particle composition 

Particles collected by impaction were individually sized and (where 
possible) morphologically identified. For the marine boundary layer 
samples the majority of particles examined could be clearly identified as 
having a marine (sea-salt) origin although occasional acidic particles and 
carbon clusters were observed. Fig. 20 shows a typical field of particles 
collected near Cowley Beach at about 1.7km and photographed using 
transmission electron-microscopy. The prevalence of cubic salt particles 
with their characteristic annular ring is clear. Another transmission 
photograph from the Cowley Beach area showing a rarer carbon cluster 
particle collected at about 900m in a maritime air stream is shown in 
Fig. 21. Approximately 12% of particles in the (marine) boundary layer were 
irregular and could not be identified morphologically. Free-troposphere 
particles exhibited a completely different character as shown in Fig. 22. 
The majority (76%) had the typical sulfuric acid or partially amiaoniated 
sulfate morphology of a central cap-shaped particle or sulfate bar 
surrounded by rings of droplets. Again a fraction of the particles, but in 
this case 24%, could not be identified in this way. For the purposes of 
calibration of the ASASP all particles in the free-troposphere were 
assumed to be ammonium bi-sulfate and for the boundary layer, particles 
with radii less than 0.3pm were assumed to be ammonium bi-sulfate and 
larger particles hydrated sea-salt. 

Particle size distributions 

Because of the failure of the Climet PHA fewer size distributions than 
originally planned could actually be determined. Large particle 
concentrations were derived by manual sizing of impaction samples which is 
a tedious process and which to a large degree was limited by the available 
resources that could be applied to the analysis. Consequently not all of 
the collected samples have been analysed (sized) . It is fairly clear from 
the observations however that particularly in the free troposphere the 
majority of variations in the aerosol involve changes in the overall 
concentration (as shown earlier, up to two orders of magnitude at least) 
with relatively minor changes in the relative concentrations of different 
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size particles, that is in the size distribution. The approach followed 
here has been to derive a number of "model" size distributions by 
combining data obtained in broad latitude ranges and a few fairly wide 
altitude ranges. The basic regions selected were 4°S-15°S, 15°S-25°S and 
35°S-45°S and the altitude ranges 0-.76km, .76km-2.3km and 3.5km-6.7km. 
These altitude ranges principally separate free-troposphere and boundary 
layer particles with one transition range. For 2.3km-3.5km the 3.5km-6.7km 
model has been used. Some further rationalisation has been made by 
combining all the impaction data into sets ( the majority obtained in the 
Cowley Beach area) differentiated only with respect to altitude and using 
the simultaneously observed concentrations of particles with radii greater 
than 0.2pm on the ASASP to merge this data with the other size data. This 
form of merging was required because the diffusion battery and impaction 
samples were not always concurrent whereas ASASP concentrations were 
obtained simultaneously with all other determinations. 

The averaged data used to derive the Cairns region free-troposphere size 
model are shown in Fig. 23 identified for their method of measurement. 
Similar averaged data for the 0.76km-2.3km altitude range and the 0-.76km 
range are given in Figs. 24 and 25. Particles with radii less than 0.005pm 
were deleted from the set for fitting an analytic distribution function. 
Distribution functions comprising two log-normals were fitted to the data, 
expressed as log( dV/dlog R) by (non-linear) least squares methods. For 
the Cairns data the free troposphere and boundary layer distribution large 
particle mode radii were constrained to 0.3pm in order to obtain 
distributions that were reasonably physically realistic. This wasn't 
necessary for the 0.76-2. 3km distribution. The separate and combined 
distribution functions are indicated with the measured data in Figs. 23 to 
25. For the 0-0. 76km data an alternate size distribution function for 
particles with radii larger than 0.22pm was derived by combining two power 
law segments fitted between the radii of 0.22pm to 0.5pm and 0.5pm-20pm 
using least squares methods. For the power law distributions shown, the 
data from the ASASP for radii greater than 1pm was excluded because of the 
possibility of an inertial loss mechanism that the rapid roll-off of this 
data suggests. In fact the removal of this data has very little effect on 
the fitted power law distribution. Model distribution parameters are 
summarised in Table 1. 

Mean distributions for Sale for 0.76km-2.3km and 3.5km-6.7km obtained in 
the same manner as the Cairns data and similarly fitted size distributions 
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are given in Figs. 26 and 27. In both cases the large particle mode radius 

was constrained to 0.3pm. 

© 

Uncertainty in aerosol size distributions 

Aerosol particle size distributions measured from an aircraft platform are 
potentially susceptable to a large number of possible error sources. Some 
of these may be systematic, such as error in sizing particles in a light 
scattering spectrometer because of unknown refractive index or uncertainty 
in the calibration relation, or loss of particles in a sampling duct. 

There is the statistical problem of collecting or counting enough 
particles to reduce the variability to an acceptable level, a typical 
problem at large radii, and the related problem of how representive are 
area, altitude or time averages used to bring up the number of particles 
which are present only in very small quantities. If an analytic function 
is fitted to the data there is the question of how well this function 
represents the data. Clearly the relative contributions of these various 
error sources can vary from instrument to instrument and even as a 
function of size for any one instrument. 

Every effort has been taken to minimise the potential for errors in this 
work. Size spectrometers were calibrated with known size particles 
before, during and after the flights and the best known calibration 
factors have been used to convert to ambient particle sizes. As an 
estimate, particle sizes should be accurate to about 5%. Because the ASASP 
was sampling from an aircraft "isokinetic" line there should not be any 
systematic inertial sampling losses however there does appear to be an 
apparent significant loss of particles with r > 1pm in the boundary layer 
with this instrument. Inertial losses certainly should not be a problem 
for particles with radii less than about 0.5pm. An alternative explanation 
(for the apparent rapid roll-off in concentration for r > 1pm) may lie in 
the calibration relation for this instrument particularly at these large 
sizes where the theoretical and experimental agreement is not as good as 
at smaller sizes. For particles with radii less than about 0.5pm count 
statistics are not a limitation in the mean distributions although there 
may be some questions regarding the use of altitude averages, particularly 
at altitudes near the inversion where size distributions do in fact change 
quite rapidly. Flow rate was measured using a mass flow sensor and should 
have less than 5% uncertainty. 
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Estimation of the accuracy of distributions obtained from inverted 
diffusion battery penetrations where there is real noise present is a 
particularly difficult problem. The GE counter was calibrated against a 
"standard" Poliak counter and concentrations obtained with it should be 
accurate to about 10%. Previous simulations of inversions with typical 
noise^ showed that uncertainty in locating the distribution peak was about 
5%-10% and in determining the magnitude of the peak about 20% in 
individual size distributions. Averaging many distributions as has been 
done here will reduce the statistical uncertainty but since not all 
distributions represent the same air parcel there is still a question of 
what the mean distribution actually represents. 

Calculating size distributions from impaction samples is particularly 
troublesome with respect to uncertainty. It is very difficult to 
accurately gauge the accuracy of any particular model used for converting 
dry particle sizes to ambient humid pre-impaction conditions in an 
arbitrary situation. The model used for sulfate particles has been tested 
previously in the stratosphere*® and is believed to be better than to 
within 20%, possibly 10% for determination of pre-impaction radius. The 
model used for converting salt particles likewise is probably accurate to 
around 10 or 20% in deriving pre-impaction radius. Jet impaction samples 
were taken from the aircraft isokinetic line so there is some potential 
there for added inertial losses. In the free troposphere a major 
limitation is the small number of particles with radii larger than about 
0.5pm, for the impaction sample means there were no particles larger than 
lpm radius, statistical uncertainty is small for most of the impactor 
concentrations, shown for example in Fig 23, but increases fairly rapidly 
for the largest particle sizes. The third largest size has a standard 
deviation of 11% in concentration, the second largest 15% and the largest 
point at apx. 0.8pm radius 33%. 

The overall good agreement between the concentrations obtained using quite 
different sizing methods as shown for example in Figs. 23 and 24 gives 
some confidence to the accuracy of the procedures and calibrations used. 

It is clear that for this type of work a diversity of measuring methods 
and instruments is almost mandatory. 
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Backscatter calculations 


Computation of backscatter from particle size distributions serves two 
main purposes. First by determining the differential backscatter as a 
function of particle radius the "active scattering” size range of the 
particle population for the wavelength of interest is readily determined. 
Second, the calculated backscatter integrated over all particle sizes 
provides a means qf comparing experimentally determined aerosol 
concentration values and experimentally determined backscatter values. 
Backscatter was determined from the aerosol particle model size 
distributions using a Mie solution program written by S. Banks and which 
utilises algorithms given by Wiscombe 11 . For particles in the large 
particle mode refractive indices used were based on the morphologically 
identified major component. For the accumulation mode a composition of 
ammonium bi-sulfate was assumed in all cases. Three wavelengths have been 
considered, 0.532pm, .693pm and 10.6pm. No values for the refractive index 
of ammonium bi-sulfate could be located for a wavelength of 10.6pm and so 
a value of 1.98-0.06i, a value reported for ammonium sulfate^ has been 
used. 

Backscatter distributions for the Cairns free troposphere distribution 
(Fig. 23) are given in Fig. 28 for the wavelengths of 0.532pm, 0.693pm and 
10.6pm for a composition of ammonium bi-sulfate. (Particle morphology 
indicated only partially neutralised sulfate as the dominant large 
particle type). Backscatter at 10.6pm is clearly dominated by particles in 
the 0.5pm to 2pm radius range whereas for the shorter wavelengths a little 
over 60% of the calculated backscatter is due to the accumulation mode 
particles (with radii less than about 0.15pm). The remaining fraction is 
due to particles with radii up to about 1pm. 

The situation is somewhat different in the boundary layer (see Fig. 29) 
where at both short (0.532pm) and long (10.6pm) wavelengths the 
backscatter is controlled principally by particles in the 0.5pm to about 
3pm radius range, although at 0.532pm there is still about a 25% 
contribution from accumulation mode particles. For the boundary layer 
calculations a refractive index for water was used for the large 
particles, a reasonable approximation in line with the observed dry 
composition of sea-salt and the high humidities (>85% RH) which were 
observed during the measurements. 


All 




Mean profiles of calculated backscatter coefficient as a function of 
altitude are given in Fig. 30 and Fig t 31. These were derived by scaling 
the concentration of particles with radius greater than 0.25pm (ASASP 
channel 35+) from the mean concentration profiles (Figs. 17-19) with the 
backscatter per particle calculated from the model aerosol size 
distributions. In both cases the free troposphere size distribution model 
was used down to 2.3km. For the mean profiles of 10.6pm backscatter 


coefficient. Fig. 30 only the Cairns region profile below 2.3km is given 
but the effect of using two different refractive indices is shown. Profile 
(a) was derived with a refractive index for water m=1.18-0.67i in the 
large particle mode and profile (b) for a "maritime" aerosol using m=1.38- 
0.057i*^. The enhanced scattering coefficient at about 4.5km in the New 
Guinea profile results from the enhanced particle concentration due to 
widespread cloud that was discussed earlier. Backscatter profiles 
calculated for 0.532pm at Cairns and 0.693pm at Sale are given in Fig. 31. 
For these shorter wavelengths the amount of scatter per particle larger 
than 0.25pm radius was nearly independent of altitude (for the three model 
altitude ranges) so discontinuities at the boundaries of the altitude 
regions are small. For the 0-0. 76km altitude range at Sale (0.693pm) the 
scatter per particle determined for Cairns (at 0.693pm) was used with the 
power law distribution and a refractive index for water. A similar model 
was used for the 0.532pm backscatter at Cairns in that altitude range but 


a profile segment for backscatter calculated using the leg normal 
distribution, again with water is also given. Calculated mean values of 
backscatter coefficient for the different wavelengths and locations 
determined from the aerosol distribution models are summarised in Table 2. 


Conclusioi 


Model size distributions were derived from free-tropospheric particle 
measurements in the radius range 0.005pm and 1pm, and for marine boundary 
layer from particles in the radius range 0.005pm to apx. 20pm. In the free 
troposphere the majority of particles were found to have sulfate 
morphology typical of sulfuric acid or lightly ammoniated sulfate. Most of 
the marine boundary layer large particles were found to have sea-salt 
morphology. The concentration of large particles (radii greater than 
about 0.2pm) was found to decrease to about the height of the inversion 
and thereafter to the maximum altitude studied remain relatively constant. 
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In the more southern latitudes (around 38°S) particle concentrations were 
exceedingly low (about one tenth of the average) in the deep south- 
westerly air stream from over the Southern Ocean following the passage of 
a cold front. Concentrations in this region also demonstrated a very large 
range in concentrations, around two orders of magnitude. Variability at 
about 17°S was considerably less and was typically within a factor of ten. 
Large variability principally between one morning sounding and several 
daytime soundings was also seen near the equator below 3km altitude. 
Calculations of backscatter using Mie theory have shown that 10.6pm 
backscatter in the free-troposphere for the regions studied is controlled 
mainly by particles in 0.5pm to 2pm radius range. At shorter wavelengths 
(0.532pm and 0.693pm) particles in the main accumulation mode with radii 
less than about 0.15pm produced most of the backscatter with a minor 
contribution up to about 1pm radius. Mean values of 10.6pm backscatter 
calculated for the free troposphere from the aerosol size distributions 
were 1*10 -11 m -1 sr -1 for the Sale area,1.5*10 -11 m _1 sr -1 for the Cairns 
area and 1.7*10 - ^^m~^sr“^ for New Guinea. 
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Table la. Summary of fitted distribution parameters for bimodal log-normal 
distributions consisting of the sum of two equations of the form: 

dN/dln r= A ( (2n) ^ In a) -1 exp{ - [ln(r/r m )]K/ [2 (Inc) 54 ]} 


Distribution 

Cairns: 


A 1 

r ml 

CT 1 

*2 

r m2 

°2 

3.5-6 

km 

122.6 

.025 

1.84 

.040 

.3 

1.49 

.76-2.3 

km 

653.7 

.043 

1.61 

.649 

.38 

1.71 

0-.76 

Sale: 

km 

375.3 

.056 

1.51 

1.37 

.3 

2.40 

3.5-6 

km 

92.1 

.027 

1.76 

.038 

.3 

1.46 

.76-2.3 

km 

104.0 

.054 

1.54 

.083 

.3 

1.92 


Table lb. Power law segments for 0-.76km Cairns distribution, r>0.22jim. 
The equations have the form : 

dN/dlog 10 r«10 [A lofflO<r)-B] 


radius A 

0.22pm-0.55um .8281 

0.55pm-20iim -2.439 


B 

.8619 

-.1129 
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Table 2. Backscatter coefficient values calculated from aerosol particle 
size distributions. New Guinea backscatter determined with Cairns size 
distribution. 


Location 

wave- 

length 

aerosol 

type 

mode 

P n m“ 1 sr“ 1 

Cairns 

3.5-6km 

0.532pm 

10.6pm 

(nh 4 )hso 4 

(nh 4 ) 2 so 4 

small 

large 

total 

2.90*10* 9 

1.75*10' 9 

1.51*10' n 

Cairns 

. 76 — 2 . 3km 

0.532pm 

10.6pm 

(nh 4 )hso 4 

maritime 

(nh 4 ) ? so 4 

maritime 

water 

small 

large 

small 

large 

large 

4.73*10" 8 

3.86*10~ 8 

7.11*10‘ 12 

7.86*10‘ 10 

2.32*10' 9 

Cairns 

0-.76km 

0.532pm 

10.6pm 

(NH 4 )HS0 4 

water 

water 

(nh 4 ) 2 so 4 

water 

small 

large 

p.law 

small 

large 

4.95*10‘ 8 

1.66*10 -7 

2.56*10”' 

6.92*10” 12 

1.55*10” 8 

Sale 

3.5-6km 

0.693pm 

10.6pm 

(nh 4 )hso 4 

(nh 4 ) 2 so 4 

total 

total 

2.27*10” 9 

1.04*10 -11 

Sale 

.76-2. 3km 

0.693pm 

(nh 4 )hso 4 

maritime 

water 

small 

large 

large 

9.43*10” 9 

3.50*10” 9 

2.38*10“ 9 

PNG 

3.5-6km 

10.6pm 

(NHO,SO a 

total 

1.71*10 -11 


A15 




A16 


Figure 1. Map of Australia and New Guinea showing principal experiment 
sites at Sale and Cairns and flight tracks. 




(«»>n Baniinv 


w a. 
o co 
— «< 

CO 

% < 
u 

M-4 S % 
O SI >1 

CO f-l 

WHO) 

a> • > 

H O -H 
-H 4-> 
'W'O U 

o c <u 

J-i <d O. 
O. «0 
6 O 
c si u 
O v£> 

•H © ^ 
4-> • > 

(d o ^ 

«M xv O 

a u 4-> 

<D 

L> U3 -v 

c 4-> -i-i 

O -H 

O * 

CO 

• co a> 

— v a> »-i 

4-> rH *H 

<3 O »W 
O -H O 
U 4-> U 

m a 
a> <d 
tj ft ‘ 


VO «-— >■ <D 

co a a 
<r» St-H 
H CJ H 

-o u 

CNJ • 

o a 

>1 Si 

Ctf *10 

5C K • 
P.O 
0) <0 
•H "O 

* ^ c 

co <d 

• a 

Csj - St 
S*l C4 
0) 0) • 
Jj 4J O 

S3 C 
t* =3 *C 
•HOC 
U «J 


A17 




(«*M) aaminv 


8 


Figure 3. Cape Grim May 30,1986 (descent), same size limits as Fig. 2. 
These profiles and those in Fig.l show two orders of magnitude change in 
concentration in the 3.5-4km altitude region over a 1 day period. 
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Figure 6. Stream analysis for 2300Z (1100, May 29 1986 local time) at 500 
HPa. Sale is circled and the stream from the Southern Ocean passing over 
Sale is indicated by the dashed line. 
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Figure 7. Stream analysis for 2300Z (1100, May 30 1986 local time) at 500 
HPa. Sale is circled and the dashed line indicates the stream line for air 
passing over Sale. This indicates a trajectory recently passing over 
central Australia. 
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Figure 8. Stream analysis for 2300Z (1100, June 10 1986 local time) at 500 
HPa. Cairns is circled and the stream passing over Cairns is indicated by 
a dashed line. 
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Figure 9, Profiles of the concentration of particles with r>0.2pm (Climet 
channel 2+) in the vicinity of Sale-Cape Grim for May 26-June 1 1986. The 
enhancement of concentrations at about 3.5km followed the reversion to a 
continental trajectory after a period of oceanic trajectories. 
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Figure 10. Profiles of the concentration of particles with r>0.2pm (Climet 
channel 2+) for soundings in central New South Wales made during transit 
flights on June 1 and June 11 1986. These profiles are very similar to 
those from Cairns. 
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LOG (CONCENTRATION) (cm 3 ) 

Figure 11. Profiles of the concentration of particles with r>0.2pm (Climet 
channel 2+) for soundings in the Cairns-Cowley Beach area, June 2 to June 
11 1986. 
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LOG (CONCENTRATION) (cm 3 ) 

Figure 12. Profiles of the concentration of particles with r>0.2pm (Climet 
channel 2+) for soundings made in New Guinea June 5 and June 7 1986. The 
profile with the minimum concentrations below 3km was obtained near Rabaul 
at around 9 AM. The other soundings were taken around midday. 
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LOG (CONCENTRATION) (cnf 3 ) 

Figure 13. Concentration profiles for particles with r> 0.13pm (ASASP) from 
Sale-Cape Grim soundings Hay 26-June 1 1986. Corresponds with Climet 
soundings in Fig. 9. 
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Figure 14. Concentration profiles for particles with r>0.13pm (ASASP) for 
central New South Wales soundings on June 1 and June 11 1986, Corresponds 
with Climet data in Fig. 10. 
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Figure 15. Concentration profiles for particles tfith r>0.13ym (ASASP) for 
Cairns-Cowley Beach soundings June 2-June 11 1986. Corresponds with Climet 
soundings in Fig. 11. 
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Figure 16. Concentration profiles for particles with r>0.13pm (ASASP) for 
soundings in Papua New Guinea June 5-June 7 1986. Corresponds with Climet 
soundings in Fig. 12. 
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Figure 18. Cairns-Cowley Beach June 2-June 11 1986. Mean concentration 
profiles, same size limits as Fig. 17. 
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Figure 19. Papua Mew Guinea June 5-June 7 1986. Mean concentration 
profiles, sane size limits as Fig. 17. 
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Figure 20. Transmission microscopy image of particle field from collection 
at 1.7km near Cowley Beach in maritime air stream. Sea salt particles are 
rectangular with low debris rings surrounding them. 
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Figure 21. Transmission microscopy image showing a carbon (?) particle 
cluster collected at 900m in a maritime airstream near Cairns. 
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Figure 22. Transmission microscopy image showing a particle field from 
collection at 6km near Cairns. Note that most particles have a sulfate 
morphology of a central cap-shaped particle (or ammonium sulfate bar) 
surrounded by droplet rings. 
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RADIUS (yum ) 

Figure 23. Mean size distribution data used to fit analytic distribution 
function for Cairns, altitudes 3.5-6km. Circles are GE diffusion battery, 
crosses ASASP, triangles impactor. The two fitted distribution modes and 
combined distribution function are also given. 




A 39 


Figure 24. Same as Fig. 23 but for Cairns 0.76km-2.3km altitudes. 
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RADIUS ( jum ) 

Figure 25. Same as Fig. 23 but for Cairns 0-.76km. As well the alternative 
power law distributions for large particles are shown. 
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Figure 26. Same as Fig. 23 but for Sale 3.5km-6.1km. 
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Figure 27. Same as Fig. 23 but for Sale 0.76km-2.3km. 
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Figure 28. Cairns 3.5km-6km backscatter coefficient distributions 
calculated from the size distribution in Fig. 24. The short dashed line is 
for a wavelength of 0.532um, the solid line 0.693ym and the dot-dashed 
line 10.6)im. 




< 1 > T 3 a> 

4-> a> c 
aJ -h 



w A) a> ^ -p 

xz 4 -> <0 

t) P « fl u 

u u «a 

0 a <o v£> j« 

& O ft • o 

*H l*i 0) O 
U4 «M « H ,Q 


11 1 1 1 I I 1 I 


•«) (y 3010/00)301 


A44 



So SI- 


i 



1 1 1 1 1 L 

hs CD ID =J- O CM 

(wm) saniiinv 


J2 

<S l 


A45 


Figure 30. Profiles of calculated 10.6um backscatter coefficient as 
function of altitude in three latitude ranges. Profile (a) used the 
refractive index of water and (b) a "maritime" value (see text) . 
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LOG $ ( (n"sr ) 

Figure 31. Profiles of calculated backscatter coefficient for 0.532pm at 
Cairns (solid line) and 0.693pm at Sale (dashed line). Profiles derived 
from model size distributions and the mean ASASP concentration profiles 
for particles with r>0.25pm. 
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Introduction 

NASA is planning a space- based Laser Atmospheric Wind Sounder (LAWS) to 
measure global wind velocity profiles. One of the most important design issues 
will be the values of aerosol backscatter coefficient on which to base the LAWS 
system design. The most uncertain region of the globe is the low latitude 
southern hemisphere where aerosol backscatter is expected to be at near minimal 
values. In order to obtain aerosol backscatter values and to assist in 
resolving this issue, NASA/CSIRO conducted an experiment in 1986 with ground 
and airborne lidars combined with an extensive set of aerosol sampling 
equipment on a CSIRO F-27 aircraft. 

A joint NASA/CSIRO flight experiment was conducted in Australia during 25 
May to 11 June 1986 to determine the aerosol backscatter at the 0.532, 0.694 
and 10.591 fan laser wavelengths. The ground-based 0.532 Jim doubled Nd: YAG 
lidar was based at Cowley Beach, Australia; the ground-based 0.694 /an ruby 
lidar was based at Sale, Australia; and the 10.591 /an coherent continuous- wave 
(CW) CO 2 lidar was mounted in the CSIRO F-27 aircraft. Aerosol equipment was 
also mounted in the aircraft which determined the aerosol number density, size 
distribution, and species. Comparison measurements of aerosol backscatter were 
made between the airborne CW CO 2 lidar and (1) the ruby lidar at Sale, 
Australia, and (2) the doubled Nd: YAG lidar at Cowley Beach, Australia. 
Several special flights were conducted by the F-27 aircraft for just the CW CO 2 
lidar and aerosol sampling measurements. 

The CW CO 2 coherent lidar was developed by NASA-MSFC and was operated by 
William D. Jones of NASA-MSFC during the flight. Jones coordinated the 
installation of the CW lidar on the F-27 aircraft, operated the lidar during 
the experiment, and guided the analysis of the data. The lidar operated 
successfully for the duration of the experiment from 25 May to 11 June 1986. 
Figure 1 shows the aircraft flight track and CW CO 2 lidar measurement path. 
The Cowley Beach location is located just south of Cairns. Measurements were 
made from 0 - 7 km altitude, from 5° S to 42° S latitude, and from 145° E to 
153° E longitude. A total of 60 flight hours were accumulated during the 2.5 
weeks of the experiment. The ''cloud-free" lidar data spans about 23 hours as 
shown in Figure 2. Measurements with the CW CO 2 lidar in the free troposphere 
gave a backscatter coefficient of 1-2 x 10 -11 m -1 sr -1 at 10.6 /im although the 
accuracy and sensitivity of the instrument were marginal at this level. 
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It is recommended that mare experimental data be taken in this region, 
that the lidar be upgraded in laser power in order to increase system 
sensitivity, and that the CW lidar data calibration methodology in both the 
volume mode and single particle node, be reviewed and improved. 



EAST LONGITUDE 

Figure 1. Flight track of CSIRO F-27 Aircraft. 
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2. Description of the CW CO 2 Lidar 

A simplified optical layout of the CW CO 2 lidar system is shown in Figure 

3. Detailed descriptions of the lidar system are presented in References 1 and 
2, and a photograph of the lidar is shown in Figure 4. The lidar uses a 10.6 
/im CO 2 laser with a typical output power of 7 W. A Mach-Zehnder interferons ter 
is used to separate the laser beam into a transmitted and a reference beam. 
The transmitted beam is directed into a 15-cm off-axis telescope which is used 
to both expand the transmitted beam and to collect the backscattered, 
Doppler-shifted laser radiation. Focus adjustment of the telescope is 
possible, varying from 10 m to infinity. A focal range of 10 m was primarily 
used for this measurement program, except at low altitudes where the laser beam 
was collimated. 

The expanded beam from the telescope is directed into the atmosphere by 
using two flat mirrors. The second mirror is able to oscillate about its 
vertical axis, thus changing the angle between the lidar line-of-sight and the 
aircraft flight direction. By oscillating between two angles during flight, 
the aerosol signal is shifted into and out of the pass-band of the processing 
electronics, due to the different Doppler frequency shifts. The processor 
bandpass filter can be tuned to allow a broad range of aircraft speeds. For 
this measurement program, the two angles were +17.44° (where 0° would be 
perpendicular to the aircraft 
position, and +14. 12° for the 


velocity vector) for the data (or signal+noise) 
noise (out of pass— band) position. 
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Figure 3. Optical Layout of the CW CO 2 Lidar System 
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Figure 4. Focused Coherent Lidar System 
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3. Theory of CW Lidar Backscatter Measurements 
3. 1 Volume Mode 

When a focused CW lidar is detecting backscatter from a target of 
"infinite" extent (e.g. aerosol particles in the atmosphere) , the theoretical 
signal- to-noise ratio (SNR) is given by: 2,3 



where VM stands for volume mode, Bp [W] is the transmitted power, h = 6.6262 
10” 3 4 [Js] = Planck's constant, V is the laser frequency [Hz], B is the 
electronic bandwidth [Hz], 7 is the optics and detector efficiency, P is the 

atmospheric aerosol backscatter coefficient [m _1 sr" 1 ], A = c/V is the laser 
wavelength [m], c = 2.9979 10 8 [m s” 1 ] is the velocity of light, f is the focal 
range of the telescope [m], and R is the e -2 (13.5%) intensity radius of the 
transmitted Gaussian beam [m] . In Eq. (1), the Gaussian beam is assumed 
untruncated and monochromatic in a coaxial system, atmospheric extinction and 
refractive turbulence effects are neglected, the local oscillator (LO) beam for 
heterodyne detection is assumed to be a Gaussian plane wave, square- law 
detection and shot- noise- limited detection are assumed, and other sources of 
received power such as molecular backscatter are neglected. Dickson 4 has shown 
that the effects of beam truncation are only negligible if the minimum physical 
lens or mirror radius is at least 2R (i. e. e~ 8 □ 0.034%). The factor 2 in the 
denominator of Eq. (1) is due to an assumed photoconductive detector. It 
should be removed for a photovoltaic detector, which is the type employed in 
the MSEC CW lidar. 

For a given set of system parameters, Eq. (1) may be written 

SNR^ = K 0 , (2) 

where K [m sr] is a system constant. If K can be determined, then P can be 
found once SNR^ is determined from measurements. Since the MSFC CW lidar 
alternately measures signal plus noise (S+N) and noise (N) through its angular 
Dicke switching, SNR^ is computed as 

S *VM = - !• < 3 > 
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A second technique for determining P employs a hard calibration target. 1 
A plane rotating target with a known bidirectional ref lectance-dis tr ibut ion 
function (BRDF) f [sr -1 ] is placed at the beam's focus. The measured SNR is 
given by 

SNR^ = G f r , (4) 

where HT stands for hard target, G [sr] is a system constant, and where £ 
depends on the illuminating and detected polarization of light and is assumed 
constant over the illumination and detection solid angles. The target is then 
removed and the aerosol backscatter SNR is given by 

SNR^ = G AL P, (5) 

where AL [m] is defined as the range interval from which 50% of the detected 
heterodyne signal is returned: 1 


AL = 



Eqs. (4) and (5) are then used to solve for backscatter: 


(6) 


SNR™, 
/? YU 

^VM " SNRjjj, 



m 

\ • / 


The advantage of the second technique is that the system constant G cancels and 
therefore many difficult to measure parameters of the system (especially fj) are 
not needed, as they are when using the first technique. Disadvantages of the 
second technique include the wide dynamic range of signals from aerosol 
particles and hard targets , and the difficulty of accurately knowing f^. s The 
second technique was used in deriving the volume mode backscatter coefficient 
for all the Australian data. 


as 


Lawrence et al 6 give without derivation AL for untruncated Gaussian beams 


0.64 ALf 
R 2 




( 8 ) 



where L ~ f in the strong focusing regime, and Post et al 7 increased AL by a 
factor of 2.2 or 


AL = 


LAM1 




(9) 


due to tbs assumption that the primary mirror truncates the Gaussian beam at 
the e~ 2 intensity points. Eq. (6) is clearly bounded by Eqs. (8) and (9). 
However, both Ref. 6 and Ref. 7 define AL as the distance between the points 
where the SNR drops to one half of its maximum value at the focus. This is not 
necessarily the same as the definition in Eq. (5). In addition, the definition 
in Eq. (5) appears to include an intrinsic factor of 2 error for SNR^ since 
only 50% of the returned energy is included. Using intermediate results in 
Ref. 6, we have derived 8 


AL = 


(LA9M1 

R 2 


( 10 ) 


for untruncated Gaussian beams, a tight focus (ffR 2 >>Af) , 
definition of References 6 and 7. (For the MSFC CW lidar, 
1. ) Eq. (10) should be compared to Eq. (8) , and is smaller 


and for the AL 
ffR 2 /(Af) * 167 >> 
by a factor of 


1.31. 


However, the definition of AL by Jones et al 1 in Eq. (5) is more 
fundamental to the calibration of backscatter. We may start with Eq. (26) of 
Ref. 3 which gives the SNR of a CW lidar focused at range f and detecting 
aerosol backscatter starting at range Li and ending at range Lz- 



Eq. (11) reduces to Eq. (1) when Li = 0 and L 2 - ®, i.e. an infinite target. 
We now ask the question: "For what path length AL - 2e centered at z = f will 
a fraction a of the total signal of Eq. (1) be obtained?" The equation to be 
solved for € is: 


BIO 



tan " 1 
- tan ” 1 



( 12 ) 


where Ckot: 1 and £>0. As it stands, this equation must be solved numerically 
for specific cases. If we again assume a tight focus, ffR 2 >>Af, and also that 
ffR 2 >>A(f+e) , then Eq. (12) becomes 


2 ^ [ifH = « B- + « n ' 1 [tF - ]] • t 13 ' 

We now assume that (ffR 2 e)/(Af 2 )< 1. Since the second tan” 1 argument is much 
greater than 1, two different series expansions of tan” 1 must be employed. 
Keeping only the lower power terms of quantities less than 1, we find 8 


AL = 2€ « "" '■■ I 2 . (14) 

Since the length AL only includes a fraction a of the return signal, it should 
be multiplied by 1/a to represent 100% of the returned signal. This cancels oc 

-i — TT»__ / 1 A \ 3 Jt _ 1 

XU EA^. ctnu yxtsxua 


, (15) 

which is identical to Eqs. (5) and (6). Therefore Eqs. (5) and (6) are 
correct, but AL is adjusted to account for 100%, not 50%, of the detected 
signal, as it should be. There is no intrinsic factor of 2 error. Returning 
to our two assumed inequalities for c, and using Eq. (14), we find that they 
require «<2ffRV(^ 2 f 2 )®17724 anc * °z2/ir = 0.64, which is quite reasonable. 

3 . 2 Single Particle Mode 

A third technique for determining ft is called the single particle 
mode . 1 - 2 • 9 ~ i 2 The lidar system is focused to a short range (tight focus) and 
the detected signal is examined for the signature of individual aerosol 
particles passing through the transmitted beam. If the lidar is airborne and 
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its beam is perpendicular to the aircraft velocity, then the particles nay be 
assumed to pass through the beam perpendicular to the optical axis. Since the 
lidar beam has a Gaussian intensity profile, the aerosol particle signature is 
expected to have a Gaussian shape. Each particle's signature will depend on 
its backscatter cross-section, c [m 2 /particle] , and on its path of penetration 
through the beam, i.e. its position both along the optical axis and transto the 
optical axis. A measurement of the peak of a particle's signature does not 
alone yield the backscatter cross-section. 

The number of peak signals of value S is given by 


N(S) = L 



A(S,ff) 


N(c) dff , 


(16) 


where L [m] is the length of the flight path, is the minimum cross-section 
which can yield the peak signal value S, A(S,c) [m 2 ] is the projected area, on 
a plane perpendicular to the particle trajectories, of the beam volume which 
yields a peak signal value of S for a particle cross-section of c, and N(<») 
[ particles /m 3 -m 2 ] is the number of particles with cross-section c per unit 
volume per unit increment. The lidar signal processor determines the peak 
height of each pulse which meets shape and width criteria and forms a histogram 
of counts vs peak height S. A large number of particles are typically used for 
each measurement, so that a uniform spatial distribution of particle 
trajectories utilizing all particle sizes can be assumed. Typically, this 
number is 10 4 . If L and A(S,Gj are known, then N (c) may be found iteratively 
until a calculated histogram N(S) matches the measured histogram. Values for 
A(S,cj are determined by spatial mapping of the system's sensitive volume. 13 
To retrieve N(cj, the following algorithm is followed: 


1. A physically reasonable model for the aerosol cross-section density 
N(Gj is postulated. 

2. The parameters of the model are postulated. 

3. For the chosen model and its parameters, the expected histogram N(S) 
is computed using Eq. (16) and the measured system sensitivity 
profiles. 

4. The calculated histogram is compared to the measured histogram in a 
least mean square sense. 
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5. 


The parameters of the cross-section density are adjusted. 


6. Steps 3-5 are repeated until a satisfactory match of histograms is 
achieved. 

The single particle node (SPM) backscatter coefficient, [m -1 sr -1 

particle -1 ] , is then computed as 


^m = -«r d<r - <17 > 

The factor of 4 ff [sr] is apparently introduced by assuming that the particle 
scatter is isotropic. 


4. Data Analysis Procedures 

Coherent Technologies, Inc. (CTI) received four magnetic tapes from W. D. 
Jones of NASA-M3FC which contained: 1) the FORTRAN source code of the MSFC 
program for analyzing the CW lidar data, 2) a calibration data file used by the 
FORTRAN program in calculating the volume mode backscatter, and 3) the flight 
data obtained during the Australian measurements. (Problems were encountered 
with the calibration data file, but a second copy sent via modem worked fine. ) 
Much of the FORTRAN code had to be adapted and/or rewritten to allow it to run 
on our Digital Equipment Corporation MicroVax II computer. All flight data 
files were combined into a single flight database. This file was written as an 
unformatted, direct access file for rapid access of data records based on date 
and time. Comparison runs were done at MSFC and CTI on selected data, showing 
near identical results, and confirming the adaptation of the code to the 
MicroVax II computer. 


The FORTRAN program was found to process the alternating data records (S+N) and 
noise records (N) in pairs. The data were collected with a 0.5 Hz switching 
rate, so that nominally 1 s data records are interleaved with 1 s noise records 
(i. e. a 2 s period). (Apparently the software does not allow a new record to 
start if the operator is in the process of entering a comment from the 
terminal. This potential bug should be investigated. ) The processing of each 
record produced a single integer for volume node backscatter V [counts] (32 
bits) and an array (histogram) of integers [counts] representing the number 
of single particle events which had a peak height in bin i of a 256-bin, 16-bit 
(0 to 65535) pulse height analyzer. The histogram included only particle 
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events which passed certain shape and width criteria. This testing is designed 
to eliminate, for example, multiple particle events and noise events. The 
number of counts in each of the bins was then summed for the data record (D) 
and for the noise record (N). The difference in these sums for record pair j 
was 

ISOM. = Ss D . - . (18) 

ISUM. was summed for each record pair j and the sum after k record pairs was 
J 

compared to an input single particle net count threshold NPART: 


k 

ISUM, = E ISUM . > NPART ? 
k j=i J 


(19) 


Typical values for NPART appeared from the code to be 1000 or 10,000. Only 
when Eq. (19) was true did the signal processor calculate and Pc^- 


From this description, it is seen that the flight length and measurement 
time (spatial and temporal resolution) for calculating Py^ and Pc^ were 
entirely determined by the threshold on the SPM counts. Running the program 
produced highly variable spatial and temporal integration intervals for the P 
calculations. Furthermore, the processing often produced negative values of 
Py j^, while was always positive. This occured in low signal regions where 
S+N and N values were nearly equal. The FORTRAN program discarded these 
negative Py^ values, which presumably biased the overall volume mode 
backs catter results to higher values. 


An example of the output values for 25 May 1986 is shown in Figure 5. 
Note that the tine interval between calculated values of P is highly variable. 


Also, several calculations of Pc^ do not have 


a corresponding ^VM plotted. 


These represent occurrences of negative values of P, 


VM‘ 


A perplexing result seen in Figure 5 is that the values of /? QPM were 


typically larger than the values of Py^ by factors of 40-90. This 


SPM 
was quite 


consistent throughout the data. Earlier comparisons 2 also showed ^SPM > ^VM’ 
but with only a factor of 4 difference. No explanation lad been found. 
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In an effort to resolve the discrepancy and to improve the data analysis 
procedures and outputs, several changes were made to the FORTRAN code by CTI 
personnel: 


1. The single particle net count threshold, NPART, was implemented as an 
input variable. 

2. A negative particle net count cutoff was implemented. Occasionally 
the particle count difference between a data (S+N) record and its 
associated noise (N) record shows a large negative excursion. A 
value of ISUM^ below -1000 causes the program to reset ISUM^ to 0. 


3. A cutoff for the elapsed tine between records was implemented. If 
the time difference between two data records is greater than 60 s, 
the value of ISUM^ is reset to 0. This prevented data that occurs 

before and after a period when the lidar is off from being combined 
into an estimate of ft. 


4. 


A cloud filter option was implemented. A flag was written to a 
data-noise record pair if the measurement was taken in a cloud. This 
was determined manually by comparing the single particle bin counts 
with the flight log from Dr. J. L. Gras. When anomalously high 
particle counts (especially in the high peak signal bins) occurred 
concurrently with reports of cloud occurrence in the flight log, the 
flag was set true. (In particular, bins 190, 192, and 256 often had 
high counts and were closely monitored. Bin 256 would contain the 
A/D output of 255 = 1111 111 IB and represents all high peak signal 
values. Similarly, bin 192 would contain 1011 1111B and bin 190 
would contain 1011 110 IB. Perhaps bits 2 and 7 were often misread as 
0 when they were actually 1. ) During execution of the program, the 


UlCLy UiJUOC »w J_V-/ i-ttX 


±A=IJL_1..3 . 


Five new data fields were also written to each record: 

1. Cloud flag. This was described above in (4). 

2. Altitude. Altitude as a function of time was computed for each day 

from data supplied to CTI by Dr. Gras. An altitude value was 

interpolated for each record by using 

ALT = 146.2 [1 - (SP/QNH) 0 - 1802 ], (20) 

where ALT is altitude [kft], SP is static pressure [mbar], and QNH is 
sea level pressure [mbar]. Altitude data were entered directly from 
the flight log for 4 June 1986 from 1105 - 1218 hours due to a gap in 

Dr. Gras' data. (1 kft = 304.8 m; 1 atm = 1013 mbar. ) 

3. Latitude. Values as a function of time for each day were hand 
entered from the flight log supplied by W. D. Jones. 

4. Longitude. Same as latitude. 

5. TAS. True airspeed values were calculated for each record from data 
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supplied by Dr. Gras using 

TAS =. [2170 (273.15 + T) DP/SP]) 0 - 5 , (21) 

where TAS is true airspeed [kt], T is corrected reverse flow 
temperature [K], and DP is differential pressure [mbar]. (1 knot = 
6080 ft/hr = 1. 151 mi/hr = 0.514 m/s. ) The TAS had been permanently 
set at 170 knots in the code. 

As an example of the usefulness of these changes, Figure 6 again shows the 

25 May 1986 plot of Ay^ vs time. The altitude of the CSIRO F-27 aircraft is 

also plotted vs tine. "Cloud" record pairs were excluded, and a value of 1000 

was used for NPART. (Recall that NPART still affects the timing of the A^ 

calculations even if A spM is not desired. ) The horizontal line segments near 

log (beta) = -5.5 represent the time intervals that the CW lidar was taking 

measurements . The calculated values of Ay^j are plotted as horizontal line 

segments of varying length. The line segments indicate the time required to 

make that measurement of ^VM' Note that there are several times when the lidar 

is on, but that no value for ^VM is plotted. This can be due to 1) negative 

values of A™, 2) ISUM was reset to zero due to reaching the negative threshold 
VM 

or elapsed time threshold, or 3) ISOM never reached the value of NPART due to 
very clean air. The last possibility is the most probable since the trend of 
decreasing Ay^ with increasing altitude is evident. 

Although these changes to the program were very useful in reducing and 
plotting the data, the large differences between Ay^ and ^SPM still remained. 
Several investigations such as varying the value of NPART were tried with no 
improvement. After several discussions with Bill Jones, we decided to 
concentrate solely on the volume mode backscatter coefficient. Plots were made 
of ^VM and aircraft altitude vs tine for the remaining 14 days of measurements 
(in addition to 25 May 1986 in Figure 6) and these are given in the next 
section. 

At this stage, negative values of Ay^ still occurred and simply were not 
plotted. We investigated whether the data (S+N) and noise (N) records could 
get out of sync with the processing program, or be reversed in order, but no 
consistent pattern in the negative Ay^ records supported this hypothesis. When 
cumulative probability plots were desired, as will be shown in the next 
section, it was not clear whether to include the effect of the negative ^VM 
values (which shifted positive Ay^ values to higher cumulative probability 
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values) or to exclude them. We calculated the standard deviation <T^ using all 
the values and then tried excluding P^ values (outliers) that were smaller 

than - k Gp , but it was hard to justify a particular value of A. 

The negative values of ^VM appeared to indicate a problem with 

establishing the instrument zero for these flights. In addition, the 

instrumental zero offset appeared to be time-varying. This is plausible due to 
the Dicke switching technique employed to zero the CW lidar. This technique 
assumes that angle tuning of the laser beam will result in either all of the 
return signal energy or none of it being passed by the detection bandpass 
filter. However, previous flights of the CW lidar took place on board the now 
destroyed NASA CV-990 aircraft. The speed of the CSIRO F-27 aircraft was much 
less than the NASA CV-990 (170 vs 390 knots), and the lidar system position in 
the F-27 prevented a larger Dicke switching angle deviation (than in the 
CV-990) from being used. These two effects combined to cause a much closer 
spacing of the return signal spectrum positions. (The angles in the CV-990 
were 7.98° and 4.65° for a differential Doppler shift of 2. 18 MHz. In the 
F-27, the angles, given in Section 2, produced a differential Doppler shift of 
920 kHz. The processor bandwidth was 800 kHz wide. ) In addition, there was no 
automatic transmission of true airspeed (TAS) data to the lidar in the F-27, as 
there was in the CV-990, forcing the operator to manually tune the center 
frequency of the bandpass filter. (The center frequency is tunable from 2.5 to 
7 MHz. ) These effects all increased the probability that the aerosol 
backs catter swas not fully transmitted by the bandpass filter in the "on" or 
"data" position, and that it was not fully blocked in the "off" or "noise" 
position. A time- varying zero offset would also be likely. 

In order to investigate this hypothesis, a set of time intervals 
representing constant altitude flight of the F-27 was selected from the flight 
database. Rather than accumulate data- noise record pairs until ISUML c NPART, 
values of P^ were calculated for each data- noise record pair. Therefore a 
varying number of ^VM calculations occurred in each time interval. For each 

time interval a mean value, was calculated. In order to examine only low 

backs catter flight segments, the value was rejected if the value of ISUM^, 
with k records in that time interval, was larger than 100. Both positive and 

negative values of ^VM were accepted. (The mean altitude value for each time 
interval was also computed for later use. ) Figure 7 shows a histogram of the 
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calc ula ted values of Note the negative mean value of these calculated 

values, and that their standard deviation equals 2.28 10 -1 1 m -1 sr _1 . Using the 

same set of time intervals, the whole process was repeated, except that data 

records were skipped, and ^VM was calculated from noise- noise record pairs. 

The resulting histogram is shown in Figure 8, using the same time intervals 

that passed the ISUM^ test for the data- noise record pairs. The noise- noise 

standard deviation was 3. 14 10" 11 m -1 sr~ l . Since approximately twice as many 

values of jflL, were calculated in the data-noise case as in the noise-noise case 
VM 

for each tine interval, we should scale the noise- noise standard deviation by 


<r = 3.14 10 _11 //2“ = 2.22 10" 11 m^sr" 1 , (22) 
eq 

which is very close to 2.28 10 _11 found from the data-noise pairs. It appears 
for this test case of tine intervals, and for ISUM^ < 100, that the backscatter 
signal is completely overwhelmed by the noise (actually the tine- varying 
instrument zero). This is supported by the positive mean value in Figure 8. 

Since the standard deviation of the calculated values of P m in the 
noise- noise case, adjusted for the greater number of data-noise pairs, is 
approximately 2 10 -11 m _1 sr _1 , this value should be used for the overall 
instrument sensitivity during the Australian flight program. It is conceivable 
that the instrument zero fluctuations were smaller (or larger) during specific 
portions of the flight program. It is highly recommended that this limiting 
factor on instrument sensitivity be investigated and improved for future 
flights. 
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Figure 5. Plot of single particle mode and volume mode backscatter coefficients. 
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F igure 7. Histogram of calculated values using data— noise record pairs. 
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Figure 8. Histogram of calculated values using noise-noise record pairs 






5. Measurement Results 

As discussed in the last section, plots identical to Figure 6 were made 
for the other days of the flight program. These are shown in Figures 9-22. 
Recall that "cloud" record- pairs are excluded, NPART = 1000, negative values of 
are thrown away, the length of each P^ horizontal line segment shows over 
what period of time it is calculated, and the horizontal line segments near 
log (beta) = -5.5 show when the CW lidar was taking measurements. Figure 23 
gives the geometric mean backscatter coefficient vs altitude for the entire 
flight program. 


O’ = ± 2.22 10 -11 m^sr' 1 

eq 


Figure 24 shows all of the calculated "cloud-free" backscatter values from 
the entire measurement program plotted vs altitude. Error bars equal to ± 1 

are plotted on each point. The right-hand error bar 
can be seen for many points that are not plotted, which results if P^ < 10” 1S , 
including the negative values of P^. (The length of the constant value error 
bars depends on the absolute position on the log P plot. ) The values of P^ 
that are clearly above the instrument sensitivity occur below 2 km and also 
near 4.5 km. Figures 25 and 26 show the data from Figure 24 divided into 
0-25 S and 25-45 S latitude bands. Figure 27 shows the cumulative probability 
plot of "cloud-free" P^ for the whole measurement program, but only for 
altitudes between 3 and 7 km. The median value of P xnj is 3 10 -11 m -1 sr -1 . 
Finally, Figure 28 shows the cumulative probability divided into four latitude 


Jis- Tra Inac ft 


IlJ3*JLJLCU.i VQ1UCD WX. 


VM 


differ by about an order of magnitude, with 


the largest value and the greatest variability (slope) occurring in the 0-15 S 
latitude band. 


6. Conclusions 

The M5FC CW CO 2 lidar was adapted for the CSIRO F-27 airplane, shipped to 
Australia, and used successfully to measure aerosol backscatter at ^ = 10.6 fm 
during a 2. 5- week measurement program in May- June 1986. The lidar data and 
data analysis FORTRAN program were transferred to Coherent Technologies, Inc. 
CTI personnel implemented several improvements to the code for reducing and 
displaying the lidar data. An instrument zero offset problem was identified 
that was tine varying and that was likely aggravated due to several aspects of 
the F-27 aircraft. Several suggestions for improving instrument calibration 
and sensitivity were made. The Australian backscatter data were presented as a 
function of time, altitude, latitude, and statistical parameters. 
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Figure 12. Volume mode backscatter and altitude vs time with clouds excluded. 
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Figure 14. Volume mode backscatter and altitude vs time with clouds excluded 
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Figure 15. Volume mode backscatter and altitude vs time with clouds excluded. 
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Figure 16. Volume mode backscatter and altitude vs time with clouds excluded. 
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Figure 17. Volume mode backscatter and altitude vs time with clouds excluded 
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Figure 20. Vo 




VJL38 901 


B36 



11 JUNE 1986 


(H»l aarnirw 

h* o in ^ n oj r> csmb 



visa son 


B37 


Figure 22. Volume mode backscatter and altitude vs time with clouds excluded 
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Figure 25. Mean backscatter coefficient vs altitude with clouds excluded 
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Figure 24. Backscatter coefficient vs altitude with - 1 (Terror bars. 
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Figure 25. Same as Fig. 24 with only 0-25 S latitude. 
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Figure 26. Same as Fig. 24 with only 25 - 45 S latitude. 
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Figure 28. Same as Fig. 27 except divided into four latitude bands. 
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Introduction 


During the period May 25-31 , 1986, CSIRO DAR operated a 
groundbased ruby lidar at Sale, SE Victoria, as part of the CSIRO/NASA 
WINDSAT experiment. The measurements were aimed at determining the c level 
of backscatter in the free troposphere above the atmospheric boundary 
layer. Accompanying measurements from the CSIRO F-27 aircraft of 
backscatter coefficients obtained in situ with a CW CO^ lidar at 10.6 ym 
and with backscatter coefficients calculated from aerosol size distribu- 
tions measured in situ are described in accompanying appendices. 

This article presents the data obtained with the DAR pulsed ruby 
lidar system together with the data reduction and analysis. 

Comparisons between the ruby and CO^ lidar backscatter 
coefficients and the calculated values from the aerosol size distributions 
are given in the main part of the report. 


2. The CSIRO Ruby Lidar System 


The CSIRO Division of Atmospheric Research pulsed ruby lidar is 
housed in a transportable caravan. The lidar consists of a Q-switched ruby 
laser mounted on the side of a 14 inch (35.6 cm) Schmidt-Cassegrain 
telescope. The system is not coaxial and maximum overlap between the laser 
beam and telescope field of view is reached at a range of approximately 
500 m. The detector is an RCA 7265 photomultiplier with an S20 
photocathode. Various combinations of neutral density filter, aperture and 
polarizer may be rotated into the receiver optics for use in specific 
applications. When the return signal is small, as it is, for instance, 
when the laser pulse is backscattered from high altitudes, the photo- 
multiplier output is amplified and electrically low-pass filtered in order 
to reduce any noise. It is then digitised by a Biomation 8100 transient 
digitiser and displayed on an oscilloscope screen. The signal returned 
from each lidar shot is saved on magnetic tape along with the relevant 
system parameters for use in the data analysis. The entire data 
acquisition process is controlled by a Hewlett-Packard HP2100 mini- 
computer and is initiated by a fire pulse from the laser energy monitor. 
Some relevant system parameters are presented in Table 1. 


Description of Data Collected 

During the period of the experiment daily measurements were made 
of the aerosol backscatter coefficient. Data were recorded twice a day for 
five days, usually in the early afternoon and evening, and when conditions 
were cloudless. A measurement was made in the morning only on the last day 
of the study. A summary of all the data taken is presented in Table 2. 

In order to collect data up to an altitude of 19 km a sequence of 
several separate series of lidar shots was recorded, each with different 
sampling intervals and trigger delays. Thus, for instance, a 10 ns 
sampling interval with zero delay was used to collect data at altitudes 
from 0 to 3 km. Similarly, the data in the range 2.1 to 8.25 km was 
measured with a sampling interval of 20 ns and a trigger delay of 13.6 ps. 
For the low altitude data no amplification of the photomultiplier output 
was necessary. For those shots which included ranges above 8.5 km, 
however, an amplifier gain of 93 was used. Unfortunately, the amplifier 
developed excessive noise during the experiment and it was necessary in 
some cases to use the low-pass filter with a cut-off frequency of 1 MHz in 
order to reduce this. Even after having done this, however, it was not 
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possible to analyse a large amount of the high altitude data due to the 
difficulty in establishing a reliable sky signal value for these shots. 

The laser energy was typically between 0.5 J and 1 J throughout 
the experiment. For the altitudes below 8.5 km up to 100 shots were 
recorded during each run. For altitudes higher than this up to 200 shots 
were recorded. Usually, one shot was recorded every 5 seconds. 

Examples of some lidar returns which were obtained during the 
experiment are shown in figures 1(a) -1(c). On the final two days of the 
experiment cloud cover often prevented any aerosol data above 1 km and 
little useful information was obtained on those days. 


Data Analysis 


During daytime hours an appreciable sky radiance D.C. signal was 
evident in the data and was measured as follows for each sequence of shots. 
The average of the last 25 data values (corresponding to the highest 
altitude) in each shot was obtained and plotted against the measured energy 
of the laser pulse. A straight line fit extrapolated to zero laser energy 
gave an estimate of the sky signal at that time. In some cases, however, 
the scatter in the data points was too great for a good straight line fit 
to be obtained and a second method was then used which relied on those 
occasions when the Biomation triggered early due to noise on the fire 
pulse. In these cases the recorded signal was due to the D.C. sky signal 
alone. 


The value obtained by either of the above methods was subtracted 
from the data for the remaining shots and the results were then converted 
into an output signal voltage using the known settings on the transient 
digitiser. The voltages thus obtained were corrected for electrical 
attenuation, for the range-squared decrease in the signal and for shot-to- 
shot variations in laser pulse energy. Voltage values were then smoothed 
over 20 m intervals. Finally, the results from all the shots in the 
measurement sequence were averaged and standard deviations calculated. 

Total attenuated backscatter coefficients were obtained in the 
following manner. The quantity measured by the lidar is: 


V T (r)* 


V T (r)r* 

E 


= K 


6 M (r) +f yr) 


] exp-2 f 

J ' n 


a M (r ’ )+a A (r,) 


dr' 


( 1 ) 


In this equation g(r), cr(r) and V(r) are the backscatter and extinction 
coefficients and the voltage induced in the photomultiplier due to 
scattering at the range r. The subscripts M, A and T refer to molecular, 
aerosol and total scattering respectively. V (r)* is the range-corrected 
voltage normalised to unit E. Pulse energy K is a system constant which 
includes such factors as the photomultiplier sensitivity, the impedance of 
the photomultiplier load, and the optical losses introduced by the receiver 
optics. In this equation it has been assumed that any contribution to the 
extinction due to molecular absorption is negligible. 
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The second measured quantity is obtained from the aircraft 
observations of the temperature and pressure profiles. These are used to 
calculate an attenuated molecular backscattering coefficient: 


att 

(r) 


e M (r) exp-2 


°M (r ' ) dr ’ 


(2) 


where 


B„ (r > * 


0.1119 c (r) 


and g (r) is calculated using the usual Rayleigh scattering equation. 
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We 

uce a calibration constant C so that at a particular range R , where 


the air is very clean, we may 

molecular and hence 8, = 0. 

0 


O' 

assume that the backscattering is totally 
We then have the relation: 




=; V ( R ) * 
C T O 


K 

C 


(R q ) exp- 2 



°M (r ' } +a A (r ' ) 


dr 


e M (R 0 )6Xp " 2 


JO 


°M (r ’ 


) dr’ 


(3) 


which allows the calibration constant C to be evaluated as: 


C 


Kexp-2 


Jo 


a A ( r’) dr’ 


(4) 


This calibration constant, which to a good approximation should depend on 
system parameters only if the aerosol loading of the atmosphere is low, may 
then by used at all the ^Jier altitudes to calculate an attenuated total 
backscatter coefficient g T (r) : 


B 


att 

T 


(r) 


V T (r)* 

~C 
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= 6 T ( r ) exp' 


- 2 j a M (r’) dr* exp-2| 


) dr* exp- 2 I a A ( r ' ) dr' 


(5) 


This is the quantity which is plotted in figure 2. We may also define an 
attenuated scattering ratio S (r) to be: 


e T (r) r r 

dr' . (6) 

and this also plotted in Figures 3. 

The practice of fitting the total backscatter $ T (r) to a 
molecular (Rayleigh) backscatter atmosphere at an altitude T where the 
atmosphere is known to be 'clean' of aerosols is standard lidar practice. 
It is only an approximate method, however, because there is always some 
a||gsol content at every altitude. We can rewrite the scattering ratio 
S att (r) 


-,att , . 
S (r) 


att, ^ 
e T (r) 

att . . 
6 m (t) 


. att , 

S (r) 


e A (r)1 


M J 

exp- 2 

J 


r r 


a A (r') dr 1 


it is found that B A (r)/g M (r) usually decreases somewhere in the troposphere 
to less than 0.1, so that C can be obtained to at least the same accuracy. 

In the present series of measurements, determination of a 
suitable normalization altitude was assisted by the availability of aerosol 
number concentrations and size distributions made in situ at a number of 
altitudes from the CSIRO F-27 aircraft (see Appendix A) . These 
measurements showed a decrease in number concentration of usually several 
orders of magnitude between 500m and about 4000 to 5000m altitude (or 
sometimes lower) . As the maximum scattering ratio was never more than 
about 5, this implied that assuming linearity between total concentration 
and backscatter coefficient, the ratio g (r)/R (r) usually fell to less 
than 0.01 at about 4000m. 1 


As the available signal Jo noise ratios in the lidar return 
degraded rapidly due to the (range) effect, (see Figures la to lc) , the 
total backscatter 6 T (r) could safely be normalised at some suitable 
altitude above 4000m without losing much meaningful information. [In fact, 
on one day, the normalization altitude was much lower. The implications in 
that case will be discussed later] . 


The normalization procedure is shown in Figure 2 where data are 
fitted to the calculated molecular scatte^j.ng atmosphere. Figure 3 
illustrates the resultant scattering ratio S (r) . It can be seen that 
the accuracy of the determination of ft A (r)/B (r) roust diminish rapidly with 
altitud^. The increase in random noise at the high altitudes due to the 
(range) ^ effects is also obvious [e.g. Figures la to lc] . 
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The derived profiles of aerosol backscatter coefficient 6 A ( r ) are 
shown in Figures 4a to 4e. A correction has been made for extinction usiny 
a value of backscatter to extinction (6/a) ratio of 6/a - 0.057 sr 
obtained from Mie calculations on the measured size distributions. As a 
first approximation, it was assumed that the backscatter coefficients found 
from the analysis of the lidar data were unattenuated and, using the 
assumed value of backscatter- to-extinction ratio, were used to calculate 
extinction coefficients at each height. From these and the measured 
molecular values the attenuation up to each height was calculated and used 
to correct the original lidar data. The process was then repeated by 
calculating new values of extinction and continued until succeeding 
iterations produced no change. The total two-way transmission to 8 km is 
about 0.65. 

In Figure 4d, the broken line indicates the result of inserting 
arbitrarily a gross error in the normalization where the data is now 
normalized to the lowest excursion in the random noise, rather to the mean 
noise. The uncertainty in 6 ( r ) due to this error in normalization is seen 
to become large above about 4 km, as expected. 

A feature on 26th May 1986 for both runs during that day [Figures 
4b and 4c] is the local sharp minimum in backscatter at an altitude of 
about 1.2 km near the top of the boundary layer inversion. The signal to 
noise ratio in the lidar return is good at these altitudes and the minimum 
appeared to be a distinct feature. It was also the point of normalization 
so that the error at the specific normalization height is again infinite 
but decreases again at the higher altitudes. 

There is some evidence for an increase in aerosol content at 
higher altitudes on 25th May 1986 and 29th June 1986, but the uncertainty 
there is high. 

There are several sources of error which have to be considered. 
The results of an arbitrary "fitting error" are shown in Figure 4d, and 
numerical estimates are given in Table 3, where the error is once again due 
to a normalization on the lowest ’swing* of the noise. Random errors due 
to the noise itself are also shown. Another error, which can be 
significant, is due to an incorrect specification of the background 1 sky 
noise*. This would have only been a significant problem in the two daytime 
runs on 26th May and 27th May. However, a similar * offset error* can occur 
on the higher gain values of the signal amplifier. This can cause the 
backscatter coefficients on various amplifier ranges to be offset from each 
other in overlapping altitude ranges. The result of offset errors is 
generally on anomalous increase or decrease with altitude depending on the 
(range) factor. Although the ratios of the gains on the various amplifier 
ranges have been normalized in overlapping height ranges, some error will 
remain. 


The uncertainty in the value of 6/a is in the region of 30%. 
However, the total two-way transmittance from 0 to 6 km is about 0.64, so 
that a 30% error causes a maximum error of about 10% in the transmittance, 
which is small compared to the other errors. 
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TABLE 1 


Ruby Lidar System 


TRANSMITTER 


Laser wavelength 

694.3 nm 

Laser energy 

0.2 - 1.5 Joule 

Laser pulse length 

60 ns 

Beam divergence 

1 mrad 

Maximum prf 

1 Hz 

RECEIVER 


Telescope diameter 

14 inch (35.6 cm) 

Receiver field of view 

10 mrad max (full angl 

Detector 

RCA 7265 PMT 

DATA RECORDER 


Biomation 8100 transient 

100 megasamples/sec 

digitiser 

100 MHz bandwidth with 


2k 8-bit words 

Computer 

HP 2100 minicomputer 
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TABLE 2 

Summary of Daily Lidar Data 


Day Time Comments 


Sun 

25/5/86 

1321 

to 

1425 

0-3 km, 1. 5-4.6 km, 2. 1-8. 3 km, 3.8-19.2 km 1 
Many mistriggers of transient recorder 

Sun 

25/5/86 

1643 

to 

1847 

0-3 km, 1.5-4. 6 km, 2. 1-8.3 km, 3.8-19.2 km^ 
Many mistriggers of transient recorder 

Mon 

26/5/86 

1217 

to 

1508 

0-3 km, 1.5-4. 6 km, 2. 1-8.3 km, 3.8-19.2 km 1 
5.1-11.3 km 1 , 9.9-31.0 km 1 
Trigger problems solved 

Mon 

26/5/86 

1914 

to 

1947 

0-3 km, 1. 5-4.6 km, 2. 1-8. 3 km, 5.1-11.3 km 1 
Night time run so sky signal close to zero 

Tues 

27/5/86 

1118 

to 

1201 

0-3 km, 1.5-4. 6 km, 2. 1-8. 3 km, 5.1-11.3 km 1 
99.0-31.0 km 

No direct sky signal measurement 

Tues 

27/5/86 

2100 

to 

2110 

0-3 km 

Night time run so sky signal close to zero 
Layer at 2.5 km altitude 

Thurs 

29/5/86 

1503 

to 

1632 

0-3 km 

Clouds above 700 metres in all of the data 

Thurs 

29/5/86 

1845 

to 

1934 

0-3 km, 1. 5-4.6 km, 2.1-8. 3 km 

Aeroplane reported very clean air above the 

inversion 

Fri 

30/5/86 

1129 

to 

1622 

0-3 km, 1.5-4. 6 km, 2. 1-8. 3 km 
Low clouds present 

Fri 

30/5/86 

1639 

to 

1733 

0-3 km 

Low clouds present 

Sat 

31/5/86 

0925 

to 

1028 

0-3 km 

Low clouds present 


Notes: 1 : The data taken for this range interval have not been 

analysed due to amplifier faults 
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TABLE 3 

Errors in 3^(r) due to incorrect normalization and due to random noise 
(background sky photon noise and signal photon noise) 


Altitude of 
normalization 

(km) 


Normalization error 

in $ (r) * 

A 



-1 

sr 


) 


Random error 
in g (r) 

, -1 - 1 , 

(m sr ) 


2 3.9 x 
4 7.1 x 
6 1.5 x 


icf 9 

8 

X 

icf 10 

icf 9 

1.4 

X 

icf 9 

10- 8 

3 

X 

ID’ 9 


* : This is due to an error of normalizing on the least 
random noise value instead of the average noise value. 


Fig. la: 

Raw lidar data taken on the 29th 
return from the 0-3 km range. 

May 

1986. 

This plot 

shows 

the 

Fig. lb: 

Raw lidar data taken on the 29th 
return from the 1.5-4. 6 km range. 

May 

1986. 

This plot 

shows 

the 

Fig. lc: 

Raw lidar data taken on the 29th 
return from the 2. 1-8. 3 km range. 

May 

1986. 

This plot 

shows 

the 

Fig. 2: 

Relative attenuated total backscatter coefficients fitted to 
molecular backscatter profile at a range of 5.5 km. Data 

the 

for 


27/5/86. 

Fig. 3: Derived scattering ratio s att (r) on 27/5/86. 

Figs.4a-4e: Profiles of aerosol backscatter coefficients (^(r) for 25/6 to 

29/6/86. The dotted line illustrates the molecular Rayleigh 
backscatter coefficient and the dashed line, (in Figure 4c) , the 
profile of g (r) obtained with normalization on the lowest 
extreme value of the random noise. (Normalized at a different 
range on 27/6 compared to the full line] . 
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Figure la 

Raw lidar data taken on the 29th May 1986. This plot shows the 
return from the 0-3 km range. 
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Figure lb 

Raw lidar data taken on the 29th May 1986. This plot shows the 
return from the 1.5-4. 6 km range. 




Figure 1c 

Raw lidar data taken on the 29th May 1986. This plot shows the 
return from the 2. 1-8. 3 km range. 



Altitude (km) 


CSIRO Ruby Lidar, Sale, Victoria 
27/5/86 1118 hrs 



Backscatter Coefficient 

(m-’sr' 1 ) 


Figure 2 

Relative attenuated total backscatter coefficients fitted to the 
molecular backscatter profile at a range of 5.5 km. Date for 
27/5/86 pi o 



© 

CSIRO Ruby Lidar, Sale, Victoria 
27/5/86 1118 hrs 



Figure 3 

n 4- 4- 

Derived scattering ratio S (r) on 27/5/86 
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CSIRO Ruby Lidar, Sale, Victoria 
25/5/86 1643 hrs 



Backscatter Coefficient 
(rrr'sr 1 ) 


Figure 4a 

Profiles of aerosol backscatter coefficients B^(r) for 25/6 to 
29/6/86. The dotted line illustrates the molecular Rayleigh 
backscatter coefficient and the dashed line, (in Figure 4c), the 
profile of (r) obtained with normalization on the lowest 
extreme value of the random noise. (Normalizied at a different 
range on 27/6 compared to the full line). 
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CSIRO Ruby Lidar, Sale, Victoria 
26/5/86 1217hrs 



Backscatter Coefficient 
(nr’sr* 1 ) 

Figure 4b 

Profiles of aerosol backscatter coefficients B^(r) for 25/6 to 
29/6/86. The dotted line illustrates the molecular Rayleigh 
backscatter coefficient and the dashed line, (in Figure 4c), the 
profile of (r) obtained with normalization on the lowest 
extreme value of the random noise. (Normalizied at a different 
range on 27/6 compared to the full line) . 
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CSIRO Ruby Lidar, Sale, Victoria 
26/5/86 1914 hrs 



Backscatter Coefficient 
(m' 1 sr' 1 ) 

Figure 4c 

Profiles of aerosol backscatter coefficients B«(r) for 25/6 to 
29/6/86. The dotted line illustrates the molecular Rayleigh 
backscatter coefficient and the dashed line, (in Figure 4c), the 
profile of B^(r) obtained with normalization on the lowest 
extreme value of the random noise. (Normalized at a different 
range on 27/6 compared to the full line). 
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CSIRO Ruby Lidar, Sale, Victoria 


27/5/86 1118 hrs 



Backscatter Coefficient 
(nr’sr 1 ) 

Fiqure 4d 

Profiles of aerosol backscatter coefficients B^(r) for 25/6 to 
29/6/86. The dotted line illustrates the molecular Rayleigh 
backscatter coefficient and the dashed line, (in Figure 4c), the 
profile of B^(r) obtained with normalization on the lowest 
extreme value of the random noise. (Normalized at a different 
range on 27/6 compared to the full line). 
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CSIRO Ruby Lidar, Sale, Victoria 

29/5/86 1845 hrs 



Backscatter Coefficient 
(nr'sr' 1 ) 

Figure 4e 

Profiles of aerosol backscatter coefficients B^(r) for 25/6 to 
29/6/86. The dotted line illustrates the molecular Rayleigh 
backscatter coefficient and the dashed line, (in Figure 4c), the 
profile of B^(r) obtained with normalization on the lowest 
extreme value of the random noise. (Normalized at a different 
range on 27/6 compared to the full line). 
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INTRODUCTION 


The lidar operated by the School of Physics and Geosciences at 
Curtin University of Technology (formerly Western Australian Institute of 
Technology) has been employed during the period 1984 to 1986 in the study 
of tropospheric marine aerosols at coastal sites in Western Australia 
(Woodman Point and Garden Island) and north Queensland (Cowley Beach 100km 
SE of Cairns) . Four field trips were conducted at the Western Australian 
sites representing temperate climates and three at the tropical site in 
Queensland. This report presents results obtained during the first 
Queensland study when comparative airborne data were obtained by the CSIRO 
Division of Atmospheric Research using the instrumented F-27 aircraft. 

2. THE CURTIN LIDAR 

The lidar uses a Quanta Ray DCR IIA frequency-doubled Nd : YAG 
laser as the transmitter and a 254mm diameter Newtonian telescope as the 
receiver. The 532nm backscattered radiation is detected using an EMI 9816B 
photomultiplier tube with an S-20 photocathode. The output signal from the 
detector is amplified using either linear or logarithmic amplifiers and 
digitised before storage on floppy discs. Digitiser operation and data 
storage are controlled via an LSI 11/23 computer. Table D1 summarises the 
lidar 1 s features. 

During operation the lidar was scanned in elevation from 0° to 
90° in a direction over the ocean. Frequently the lidar was operated at 
fixed elevation angles and 50 shot averages were recorded to improve the 
signal-to-noise ratio. Occasionally data were recorded with a polariser 
inserted in the receiver and aligned alternately parallel and perpendicular 
to the plane of polarisation of the transmitted laser pulse. 
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The lidar data were complemented by on-site measurements of 
meteorological and aerosol parameters. Temperature, relative humidity and 
wind speed and direction were recorded at ground level and profiles of 
temperature and relative humidity were obtained from radiosondes launched 
from the site during each three-hour observational period. A MRI nephelo- 
meter, a CSASP particle counter and a radon counter provided surface data 
on aerosol extinction, number, size and origin. 

3. METEOROLOGICAL CONDITIONS 

The period of observation (from 29th May to 11th June 1986) was 
characterised by SE (on shore) winds and an air temperature which varied 
between 22° and 24°C during the measurement sessions. During the earlier 
part of this period (29th May to 7th June) winds were consistently between 
7 and 11 m.sec \ relative humidities between 80% and 96% and there were 
frequent rainshowers from low cumulonimbus clouds. Conditions moderated 
during the last three days? winds dropped to between 3 and 5 m.sec 
relative humidities dropped to between 50% and 80%, clouds only formed 
after midday and there was no rain during the measurement sessions. 


4. METHODS OF DATA ANALYSIS 

The lidar data have been analysed to produce several different 
outputs. Data averaged over many shots and acquired from horizontal 
firings in regions of horizontal atmospheric homogeneity have been analysed 
to produce values of the aerosol backscatter-to-extinction ratio (B /o ) in 

A A 

the near-surface layer. Data acquired at higher elevation angles have been 

analysed to produce profiles of backscatter and extinction. In the absence 

of measured values of (S /o ) throughout the troposphere, the value 

A A 

measured at the surface was used in the analysis of these latter data. 
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4.1 The Determination of the Aerosol Backscatter-to-Extinction Ratio 

The lidar equation for the backscattered power received from a 
range r can be written as 

P(r) - K(B (r) + B (r) ) exp 
Am 

’ o 


- 


0 . (r'J 

A 


o (r’) 
m 


dr* / r 3 


(Dl) 


where B(r) and o(r) are, respectively, the volume backscatter function and 
volume extinction coefficient and K is a system constant. The subscripts A 
and m refer to the aerosol and molecular contributions to these quantities. 
Note here that the extinction terms are assumed to be due to scattering 
only; molecular and aerosol absorptions at the wavelength used are 
negligible (Elterman (1964), McClatchey et al. (1972)). 

For a lidar firing horizontally through a homogeneous atmosphere 
Dl can be rewritten ■ 

P(r) r 2 = k£ t exp[-2 o.r] , 

where 

B = B + B 

P T A p m 

is the total backscatter function and 


(D2) 


(D3) 


= o, + o 
A m 


(D4) 


is the total extinction coefficient. The bars indicate averaged values, as 
small variations about a mean are acceptable. Taking natural logarithms of 
D2 gives 
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In [P (r) r 2 ] 


(D5) 


= In [kB t 1 - 20. r . 

D5 is a linear function of r and a linear regression of In [P (r) r 2 ) 
against r yields values of o and kB^,. 

Using D3 and following a method similar to that of Waggoner 
et al. (1972) 


KB„ = kB. + kB 


m 


= K 


4tf 


A m 


(D6) 


Here the aerosol extinction and backscatter are related by the equation 


P. (*) 
A 


4tt 


(D7) 


where P (ir) is the aerosol phase function for the case of backscatter (e.g. 
A 

Deirmend j ian , 1964) , and the quantity in brackets is equivalent to the 

aerosol backscatter-to-extinction ratio. Note that this value is averaged 

over the ensemble of particles in the scattering medium. 

If D6, which is equivalent to equation 6 of Waggoner et al . , is 

normalised by dividing by 6 , then measurements made over an extended time 

m 

period or range of altitudes will all provide a linear graph of KBgp/ 8^ 

against a /8 , provided that the average backscatter-to-extinction ratio 
Am 

remains constant. D6 then becomes 


. J 


P (TO 

ace) /e = k — r- — 
T m ^ 4 tt 


o/£ + K 

A m 


(D8) 
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and linear regression will produce values of both the lidar system constant 
K and the aerosol backscatter-.to-extinction ratio. 

4.2 The Determination of Profiles of Vertical Backscatter Function 

The clean-air calibration method has been used in the analysis of 

the data. In this method, the backscatter signal at some calibration range 

r is assumed to be due entirely to scattering from air molecules, 
cal 

Because of the lidar geometry, it is not feasible to determine the trans- 
mittance from zero range to r , and it is convenient to combine the 

C3 X 

system constant and the transmittance in one term [K T 2 (0,r c _^)]. 

For ranges r < r , D1 can be rewritten 
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This can be expressed in numerical form: 
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3 3 
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Here 6 r is the digitiser range increment and P = 1.5 is the molecular 

m 

phase function for backscatter. A solution for 8 ^ is found by calculating 


first the constant term at r 


jcal 
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then using this to calculate solutions at ranges r^ — r j ca ^ “ 
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Note that the solution B. . is required in the summation on the right hand 

AD 

side of the equation. To overcome this problem this term in the summation 
initially is set to zero and a first approximation for B is obtained 
using D12. A better estimate is obtained by using this first solution in 
the summation. Iteration proceeds until the solution converges to the 
desired accuracy (Gambling and Bartusek, 1972) . 


4 -3 An Estimate of Uncertainties in the Derived Aerosol Backscatter 
Profiles 

The error (or uncertainty) in $ in D12 is found by first 
calculating the relative error in the total backscatter: 
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where AP^ arises from ^ LSB in the digitiser ADC, Ar is the uncertainty in 
the range = *5 the range increment 6r, 
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and 
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The relative error in the molecular extinction term, p, is taken as one 
percent (Lenhard, 1973) . Next A6 . is calculated as 
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The error in the aerosol extinction to be used in D14 at the next range 
increment is 




AB 


hi 
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5. RESULTS OF OBSERVATIONS AT COWLEY BEACH - MAY-JUNE 1986 

5.1 Aerosol Backscatter-to-Extinction Measurements in the Surface Layer 

Information on the aerosol backscatter-to-extinction ratio in the 
region about 2m above the water and approximately 400m to 4000m from the 
lidar was obtained from horizontal firings of the laser over the ocean. 
The lidar was situated at approximately 50m from the water’s edge. 
Observations were made several times each measurement session. 

Data have been divided into two groups for analysis. These 
correspond to the earlier rainy period of high humidity and the later 
clearer and drier period described in Section 3. As could be expected, 
data collected during the latter period show less scatter, and have been 
used for the determination of the lidar system constant K, although the 
value obtained from the first period is not significantly different. 
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Data acquired in the second period during times of horizontal 

homogeneity are presented in Figure 1. The 12 data points shown represent 

data collected at 18 different times during the 4 day period. Each is the 

average of 50 individual laser firings covering about 25 seconds giving a 

total of 900 laser firings. Note that there is remarkably little variation 

in the aerosol optical properties during this period. The results for 

P.(ir)/4ir and K are: 

A 

P A (ir)/4ir = (0.0189 t 0.0022) = (0.238 ± 0.027) / 4* sr" 1 , 

K * (3.56 ± 0.19) x 10 10 m 3 A.J -1 . 

The lidar system constant is the product of the optical transmission 
efficiency of the receiver, the effective receiving area of the mirror, the 
velocity of light and the conversion efficiency of the PMT which had been 
measured previously. The optical efficiency is the product of three 
reflections from aluminium, transmission though four glass surfaces and 
through a narrow band filter: 

i.e. K *= 


Excellent agreement exists between the directly and indirectly derived 
values of the constant. 


h S c/2 Q 
o 

(0.247) (0.0478 m 2 ) (1. 499x10® m.sec" 1 ) (2.0xl0 4 A.w” 1 ) 
3.54xl0 10 m 3 . A. J** 1 . 
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Data acquired during the earlier period are presented in 
Figure 2 which, although it shows more scatter than in Figure 1, still 
indicates rather consistent aerosol properties during the five days of 
observation. The twelve data points consist of 24 averages of 50 shots or 
a total of 1200 laser firings. Linear regression of the data gives values 
of 


P a /47t = (0.0128 ± 0.0035) = (0.161 ± 0.044) / 4* Sr 
and K = (3.59 ± 0.48) x 10 10 m 3 A.J _1 . 

Note that the system constant derived from this data set is not 
significantly different from the other values. This proves that the 
calibration of all of the components in the lidar system remained stable 
during the 14 days of the study, despite the harsh environmental conditions 
of strong winds, heavy rain and sea spray. The components considered here 
are several optical surfaces, optical filters, amplifiers, ADCs and the 
laser output energy monitor. 

The values for the aerosol phase function derived from the two 
data sets, .019 and .013, differ although the difference is only just 
significant at the one standard deviation level. An estimate of this value 
may also be obtained by using Mie scattering calculations and particle- 
sizing and nephelometer data. Total extinction at 1130 a.m. on the 

10th June was 7.3 x 10 ^ ^ which gives an aerosol extinction of 
6 x 10 ^ m . a log-normal distribution was found to fit the particle 
sizing data and when radius limits of r = 0.2 \m and r 2 = 5.0 ym were used 
with a refractive index appropriate to sea salt in Mie scattering 
calculations a value of 0.017 was found for the aerosol phase function at 
backscatter* This value compares well with the lidar derived values. 
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Finally, lidar derived values of aerosol extinction on plotted 
against nephelometer values in Figure 3. The scatter in the data points is 
the result of the instruments not sampling the same volume of air. The 
nephelometer was on the shore whereas the lidar was sampling air up to 4 km 
out to sea. The average ratio of the lidar to the nephelometer value is 
.73 with a standard deviation in this mean value of .05. As the lidar data 
were measured at 532 nm and the nephelometer at 470 nm a wavelength 
exponent of -2.6 is indicated. 

5.2 Vertical Profiles of Aerosol Backscatter Function 


Vertical profiles of aerosol backscatter function derived from 
lidar measurements on June 9th and June 10th are presented in Figure 4 and 
Figure 5 respectively. These data are studied here because airborne 
sampling of aerosols on these days allow a comparison to be made. 

Figure 4 shows enhanced aerosol scattering in the boundary layer, 
where high relative humidity has increased the particle sizes. The lidar 
profile is an average of 50 firings at around 1030 EST and the molecular 
scattering profile was obtained from a radiosonde balloon launched from the 
lidar site at 1015 EST. A weaker scattering layer is apparent between 
1400m and 2000m, also a region of high relative humidity. The aerosol 
backscatter function decreases steadily with height above this altitude 
which corresponds to that of the strong temperature inversion. Note that 
the lidar data have been normalised to the molecular backscatter signal at 
4500m and this causes the graph to tend to zero near this height. 
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Figure 5 is the average of 50 lidar firings around 1138 EST. The 
radiosonde balloon was launched at 1052. Again, enhanced scattering from 
aerosols is detected below the temperature inversion at 2100m, a region of 
high relative humidity. Above this height aerosol scattering drops 
steadily. Normalisation to molecular scattering occurs for heights above 
4000m. The aircraft particle sampling data also show the strong layer with 
weaker scattering above although, in contrast to the lidar data, they show 
roughly constant scattering with height above the inversion. The 
consistent difference between the two instruments in the measured magnitude 
of aerosol scattering is unexplained at this stage. 

Note that lidar derived aerosol backscatter function profiles are 
influenced by two factors. The first is the need to normalise data at some 
height where aerosol scattering is assumed to be below the limit of 
detection of the lidar. If, in fact, significant aerosol scattering does 
exist this region, then the lidar backscatter profiles will underestimate 
the actual value at all heights in the profile. Such effects are likely to 
be small in this data set, a fact supported by the airborne measurements. 

As discussed in section 4, a value of the aerosol scattering 
phase function in the backscatter direction (i.e. the backscatter- to- 
extinction ratio) is needed to calculate the transmittance correction from 
the backscatter measurements. If the assumed value is in error or changes 
with height then errors in the backscatter profile will result as shown in 
the analysis in section 4.3. This is not expected to be a major factor 
with data presented here. 
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Table Dl: Curtin Lidar Technical Details 


Transmitter : 

Laser 

Wavelength 
Energy per pulse 
Pulse duration 


QuantaRay DCR IIA Nd:YAG 
with frequency doubling 
532 run 
'v 0.2 J 
5 ns 


Steerable Receiver : 
Telescope 

Detection filter bandwidth 
Receiver optical efficiency 
(including filter) 
Detector 
PMT conversion 


254 mm diameter Newtonian 
1 nm 

0.247 

EMI 9816E (S-20) photomultiplier tube 
2.0 x 10 4 A.W _1 


Recording Electronics : 

Amplification 

Digitiser 

Minimum sampling interval 
Computer 


linear or logarithmic 

1024 words @ 8 bits resolution 

50 ns 

LSI 11/23+ (controls digitiser and 
recording of data onto 8" floppy discs 
Also used for data analysis) 
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